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FOREWORD 

The PLANETSYS computer progra^i (20) was originally developed 
under U.S. Army Contract DAAD05-74-C-0747 , sponsored by the B^i,ll- 
istics Research Laboratory (BRL) at Aberdeen Proving Grounds with 
Mr, W. Thompson as Technical Monitor. The PLANETSYS program sim- 
ulates the thermo -mechanical performance of a multi-stage plane- 
tary power transmission. 

This user's manual describes the use of a version of the 
PLANETSYS computer code developed under NASA-Lewis Research Center 
Contract^ with Mr. H. Coe as Technical Monitor. The revisions 
made to the program include : 


Implementation of two versions of code within the 
program- the SKF version and the NASA version. 
The differences between the two versions reside 
in the calculation of the elastohydrodynamic 
(EHD) film thickness and traction forces which 
develop between rolling element-raceway and 
rolling element-cage concentrated contacts. The 
original film thickness model CDowson-Higginson 
(13) ) and the Tallinn traction model (11) are 
used in the SKF version, while the NASA version 
uses the Loewenthal model (21) to calculate film 
thickness and the Allen model (22) to determine 
traction forces. 


Additionally, a new subroutine, FLMFAC Creplacing LRHS) , is 
used to determine the lubricant life factor as a function of A 
(film thickness/surface roughness) . The values of the lube life 
factor produced by FLMFAC adhere closely to the curve recommended 
by the ASME (23, 24) . 

A new section. Appendix F, describes the 'differences between 
the SKF and NASA methods of calculating film thickness and trac- 
tion forces, and explains the differences in executing each ver- 
sion of the code. 
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1. rNTRODUCTION 

In, a planetary power transmission system, there are usually 
three or four planet gears, equally spaced around a sun gear and 
encircled by a ring gear. In airborm^t transmissions it is common 
practice to support the planet gears 0f\ single or double row roller 
bearings such that the planet gear and bearing outer ring comprise 
a single component. The planet bearings are connected by a planet 
carrier at the bearing inner rings » Depending on the direction of 
the power flow, there are six possible kinematic inversions, which 
the planetary system may adopt. These are described in Table 1. 

The major function of the program is to compute for any of 
the six kinematic inversions, performance characteristics of a 
planet bearing. The latter may be a cylindrical or spherical 
roller bearing with an outer ring rigid or flexible. The bearing 
may contain one or two rows of rollers and its inner ring is taken 
as an elastic solid. 

The Main Program consists of the following major subprograms: 

1) Bearing Analysis - Largely based 
upon me tho ds 0 £ L iu and Chiu (1,2), 
the bearing outer ring is either 
rigid or flexible. Effects of the 
latter are reflected in fatigue 

life calculations and load distributions, 

2) Bearing Dimensional Change Analysis - 
Based on the methods of Timoshenko, (3) , 
and adapted to the bearing system by 
Crecelius, (4). 

3) Generalized Steady State and Transient 
Temperature Mapping and Heat pissipaH'on 
Analyses - based on the methods of Harris , 

(5J , Fernlund, (6), and Andreason, (7). 

In this version of the program, the user can select which of the 
following two models are employed in the calculation of EHD film 
thickness and concentrated contact traction: 

I . SKF Model 

1) An elastohydrodynamic CEHD) film thickness model that 
accounts for i) thermal heating in the contact inlet 
using a regression fit to results obtained by Cheng (8) 
and ii) lubricant film starvation using theoretical 
results derived by Chiu (9) , 
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2) A semi- empirical model for fluid traction in an EHD 
contact (10) , is combined with an asperity load 
sharing model developed by tallian (11) to yield a 
model for traction in conceatraV.ed contacts that 
reflects the state of lubrication as it varies from 
dry, through partial EHD to the full EHD regime. 

II. NASA Model 

1) In this empirical model, an elastohydrodynamic (EHD) 
film thickness is determined from the equation developed 
by Loewenthal (21). As in the SKF model, film reduction 
factors due to contact inlet thermal heating and starva- 
tion effects are applied, 

2) A model for fluid traction in an EHD concentrated 
contact was developed by Allen (22) , based on the 
fluid shear stress distribution acting over the contact. 

Both versions of PLANETSYS also make use of the following 

bearing related models ; 

1) Three components of bearing cage related friction 
are treated for both lubricated and dry operation. 

For dry operation coulomb friction is assumed. 

With lubrication, hydrodynamic' lubricant shear 
models are used as follows: 

a) Gage web-roller friction 
according to methods of 
Dowson (13) 

b) Cage rail -roller end using 
the methods for a plain 
bearing from Marks (14) 

c) Cage rail -ring land using 
the methods for a short 
journal bearing according 
to Dubois and Ocvirk (15) . 

2) A model for the effect on bearing fatigue life of _ 
the ratio of EHD plateau film thickness to composite 
surface roughness is incorporated (23, 24) . 
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Additionally, models for temperature viscosity and pressure 
viscosity variation as functions of temperature given by Walther, 
(19) and Fresco, (12) respectively, were adopted, 

The purpose of the program is to calculate the steady state 
or transient thermal performance characteristics of an airborne 
type, planetary power transmission system, comprised of up to 
three stages. Each stage consists of a unique sun, ring and us- 
ually three or four planet gears plus a carrier .. 

To save weight, the typical airborne planetary is designed 
such that the planet gear and its support bearing form a single, 
non-separable piece of hardware, as shown in Pig, la. 

It has been demonstrated in ref. (1 and 2) that the gear 
tooth loading can affect the bearing rolling element-raceway load- 
ing due to the flexibility of the planet gear. This in turn af- 
fects the heat generation rate of the planet bearing and thus the 
thermal performance of the entire system. 

The program is structured to handle up to three planetary 
stages. The power flow through each stage is assumed to be iden- 
tical. To account for the uniqueness of each stage, a spearate 
set of input data is required to describe the components of each 
stage. The total power through the stage is assumed to be 
shared equally among all planets and all planets are assumed to 
behave identically. 

Each stage may operate according to one of the six kinematic 
inversions presented in Table 1. Note that the program can handle 
each inversion and properly accounts for all planet gear and rol- 
ling element centrifugal forces which arise from the kinematics. 

The planet system thermal performance, loading, frictional 
heat generation and dimensional stability are all cross-coupled 
within the program in order to produce a complete picture of sys- 
tem performance . 

The types of planet bearings which can be analyzed are shown 
in Fig. lb. 
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TABLE I 

PLABETART INVERSION INDICES 


Inversion No. 

Power Input 

Fixed . 

Power Output 

L 

Sun 

Carrier 

Ring 

2 

Ring 

Carrier 

Sun 

3 

Sun 

Ring 

Carrier 

4 

Carrier 

Ring 

Sun 

5 

Ring 

Sun 

Carrier 

6 

Carrier 

Sun 

Ring 
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SCAB TOOTH LOAD 



OUTER BING (INTEGRAL GEAR) 
BOLLING ELEMENT 

INNER BING 

* 

i 

SHAFT 


FIGURE la PLANET GEAR BEARING SIDE VIEW 
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Single Row Spherical 



I 

ri\ 

R 




Double Row Cylindrical Single Row Cylindrical 

FIGURE lb PLANET BEARING TYPES WHICH CAN BE ANALYZED WITH 
"PLANETSYS" 

7 

SKF TECHNOLOGY SERVICES 

INDySTRiEjJNQ:; - 



oriqinai. “ 

OF POOR QUAUTT 


AT81D044 


2. PROGRAM INPUT 

2.1 Types of Input Data 

A complete set of input data Comprises data 
of 3 distinct categories* Within these 
categories, cards which convey specific 
kinds of information are referred to as 
card types. The categories are listed 
below. 

I , Title Cards 

A title card plus a second card which 
provides the program control information 
for the planet-bearing solution. 

II, Bearing Data Cards 

A set of up to twenty (20) card types 
describing one bearing in the assembly. 
All bearings in the planetary stage 
are treated identically. 

III, Thermal Data Cards 


A set of up to nine (9) card types 
to describe the thermal model of 
the assembly. 

The input data instructions are given in Appendix 
A, and are for the most part, self explanatory. The 
units used for input data are also given in Appendix 
A. They are presented in an eighty column data format, 
A description of the variables is given in the input 
instruction forms. 

2 . 2 Data Set I - Title Cards 

2.2.1 Title Card 

This card should contain the computer 
run title and any information which 
might prove useful for future identi- 
fication. The full eighty (80) 
columns are available for this purpose. 

The title will appear at the top of each 
page of Program output. 
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2.2.2 Flag Set Card 
ITEM 1 

Nun.\ber o£ reduction, stages , can be 
1 thru 3. 

ITEM 2 

Print Flag (NPRINT) , equal to zero 
is normal and will result in no 
intermediate or debug output. With 
a value of one, a low level inter- 
mediate print is obtained at the end 
of each planet- bearing iteration. 

The values of the variables, the 
inner ring displacements and the 
equation residues are printed. 

At the end of each bearing iteration, 
wherein the rolling element and 
cage equilibrium equations are solved, 
an error parameter is printed which 
has the value: 

Error Parameter » 

AXj. is the change in the variable 
X specified at iteration N. 

Xw_i is the value of the variable 
specified at the previous 
iteration. 

The Error Parameter is calculated for 
each of the bearing variables, but 
only the largest one is printed. 

Additionally, at the end of each 
Clearance Change iteration, the 
clearance change error parameter 
is printed. This error is defined; 

Error Parameter « 

Roiling Element Diamete 

Where DCL„ and DCL^,_, denote the clearance 
changes calculated‘^at the current and 
previous iterations respectively. 
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If NPRINT is set at 2 all the above 
information is printed. Additionally, 
the variable values and residue values 
are printed for each iteration of the 
rolling element equilibrium solution. 

ITEM 5 

Fit calculation flag (ITFIT), ITFIT 
controls the number of iterations allowed 
to satisfy the bearing clearance change 
iteration scheme. If ITFIT. is set to zero 
(0) or left blank, the clearance change 
portion of the program is not executed. 

If a positive number is input, the 
clearance change scheme is utilized 
with a maximum iteration limit of five 
(5). If a negative number is input, the 
scheme is used with a maximum iteration 
limit equal to the absolute value of the 
negative number. 

ITEM 4 

Main Loop iteration flag (ITMAIN) , ITMAIN 
limits the number of iterations attempted 
during the solution of the bearing equi- 
librium problem, i.e. establishing the 
equilibrium of bearing reactions and 
applied gear tooth and inertial loads. If 
ITMAIN is set to a positive integer, then 
(20) iterations are permitted. If ITMAIN 
is set to a negative number, the number of 
iterations is limited to the absolute 
value of that number. If ITMAIN is set to 
zero or left blank, only one iteration is 
performed. 

ITEM 5 

Fit Loop Accuracy (EPSFIT) , EPSFIT is the 
convergence criterion for the diametral 
clearance change portion of the analysis. 
As mentioned under item 2 above, this 
error parameter is defined to be: 
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Error Parameter ■ 

Rolling Element Diameter 

The iteration scheme is terminated when 
the error parameter is less than the input 
value of EPSBIT. If EPSFIT is left blank 
or is set to zero (0), the program default 
value of 0.0001 is used, Xn the above expres 
sion» OCL is the change in bearing clearance, 
and N refers to the iteration number, 

IT^ 6 


If Plot Flag is set at either "T*’ or 
’’.True.*’ a plot of radial deflection and load 
vs, azimuth angle will be produced at the 
end of the program output. Default is false 
(No plot). 

ITEM 7 

Units Flag (IMET) , IMET permits input 
data to be in either International or 
English units. A value of 1 indicates 
International, 0 or blank English units. 
Dimensions of specific variables are 
given in Appendix A. 

ITEM 8 

Material property Signal (IMT), a value 
of 1 allows input of material constants. 

If IMT is 0 or blank, the program assumes 
the material properties of steel apply. 

The default values are listed in sections 
2.3.14 thru 2.3.17. 

2.3 Data Set II - Bearing Data 

Most of the input instructions are self-explanatory. 
Where more explanation than given in the input data 
format instructions is required items are treated 
on an individual basis by card type and item niunber. 
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2.5.L Data Card No» 3 » Rolling Element Information 
ITEM 1 

Bearing type, columns 1-10 must be specified, 
left justified, i,e., or "C*’ in column 
This format must be followed since the Program 
recognition of bearing type, (Spherical or 
Cylindrical rolling bearing) , is derived 
from reading the ”S" or in the first column 
of this card. ' 

ITEM 8 

Roller Skew Angle, columns 71-80, is used to 
study the effect of roller skew on bearing 
generated heat and is not used in calculation 
of bearing equilibrium. Input is in 'degrees. 

2.3.2 Data Card No . 4 - Bearing Data 

ITEM 4 

Refers to the number of slices into which 
the roller raceway contact may be divided. 

A maximum value of twenty, (20) is permitted. 

A default value of two (2) is used if Item 4 
is blank or zero. Each slice is the same 
width. 

ITEM 5 

Columns 71 thru 80 contain a signal, termed 
the crown drop flag, which specifies whether 
the roller-race crown drops will be calculated, 
or read directly. If item S is blank or zero, 
the crown drops are calculated based on the 
roller-race crown radius, and flat length 
input informat ion If the crown drop flag 
is other than zero or blank the nonuniform 
separation of the roller and raceway must 
be specified at the center of each slice. 

The slice widths are identical. The non- 
uniform roller- raceway separation is input 
on data cards 7 and 7A. (see Fig. 2) 
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2.5.3 Data Card 5 - Inner Raceway Information 
ITEM 1 

Columns 1-20^ ’'Steel Designation" for the 
inner ring. Alphameric- literal description 
o£ the steel type such as "M-50" or "AISI 
52100" is input. 

ITEM 2 

"Effective Contact Length" refers to the 
longest possible length which can be ob- 
tained at a roller-race contact. Typically 
this is the roller total length less the 
corner radii. If, however, the raceway 
undercuts are exceptionally large so that 
the track width is smaller than the effective 
roller length then the track width Should 
be input. 

ITEM 5 

Columns 51 thru 60, the number input for 
item 5 is used to account for improved 
materials, multiplying the raceway fatigue 
lives as determined by Lundberg-Palmgren 
methods. Typical life factor values for 
modern steels are in the neighborhood of 

2.0 to. 6.0. In the ASME Publication 
Life Adjustment Factors for Ball and 
Roller Bearings , the tMaterial factor D 
and the Material Process Factor E should 
be used multiplicatively as inputs for 
item 5. The program computes a lubricant 
life factor based on the value of A (EHD 
plateau film thickness divided by RMS com- 
posite surface roughness). The lube life 
factor is calculated by sub-routine FLMFAC 
using a least squares fit to the curve 
recommended by ASME C23, 24) , This factor 
ranges in value from ‘0.21 for A < 0.6 to 

3.0 for A > 10.0. " 
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2.3.4 Data Card 5A - Outer Racv^way Information 

Data Card 5A, is th.e same as 5, but outer 
raceway information is given. ' 

2.3.5 Data Card 6 ^ Surface Data 

Items I ti^rough 6 define the statistical 
surface microgeometry parameters of the 
rollers and raceways . Items 1 through 
■3 require the input of center, line average 
(CLA) surface roughness. Within the 
program CLA values are converted to RMS 
by multiplying by 1.25. 

Items 4 through 6 are RMS values of the 
slopes, (degrees) of the surface asperities. 
These are measured in a traverse across 
the groove for rings, longituidinally for 
rollers. Typical values for raceway and 
rolling element surfaces are 1 to 2 degrees. 

2.3.6 Data Cards 7 § 7A - Crown Drops 

These cards are used to input the separation 
between the inner and outer raceways and the 
roller at the center of each slice along the 
roller profile with the high points of the 
roller and race in contact, i.e. with all 
clearance between roller and raceway removed. 
These cards must be omitted if item 5 of 
Data Card No. 4 is 0 or blank. 

2.3.7 Data Card 8 - Cage Information 
ITEM 7 

Assumed cage slip, input as a positive number 
less than 1, is used in the calculation of 
roller-inner raceway rolling element speeds. 
As the cage slips, sliding velocity at 
the contact is increased. Typical slip 
values range from 0,0 to 0.15. 
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2*-3*S Data Card No« 9 - Diametral Clearance 
Iteration Data 

ITEM 1 


Number of Clearance Iterations, although 
net used ia the calculations, is an input 
item. It should be set at 1., and is 
included only for future parametric studies 
of system performance as a function of 
bearing diametral clearance. 

ITEM 2 


Change in clearance for next iteration, 
set at 0., is not included in any 
calculations, (see ITEM 1). 

2,5.9 Data Card No. 10 - Planet Gear Data 


ITEM 2 


Planetary Inversion Index, refers to the 
specific kinematic condition under which 
the system operates. Valid numbers are 
one (1) thru six (6), the inversions are 
as follows: 


Inversion No. Power Input 


Fixed Power Output 


1 

2 

3 

4 

5 

6 


Sun 

Ring 

Sun 

Carrier 

Ring 

Carrier 


Carrier 

Carrier 

Ring 

Ring 

Sun 

Sun 


Ring 

Sun 

Carrier 

Sun 

Carrier 

Ring 


Note: For planetary inversions 3-6, where the 

planet carrier is not fixed, the inertial 
load of the planet gear is assumed to be 
totally taken by the planet carrier so that 
the gear tooth load remains unaffected. 
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ITEM 6 

Moment of inertia of outer ring cross section 
Cl) is calculated for input relative to the 
centre idal axis fRc) as shown in Fig. 3, The 
shaded area in Fig. 3 indicates the effective 
area used in the calculation of X. In most 
cases the neutral axis and centroidal axis are 
coincident and standard equations which relate 
I to cross section geometry can be used. The 
user may specify that the outer ring be rigid, 
(with respect to radial deflection) 'by inputting 
the moment of .inertia greater than 10^ XN^ 

(4.2 X 10® CM"). 

2.3,10 Data Card No, 11 - Misc. Information 
ITEM 1 


Power thru stage refers to the total power 
through the stage . Internally, the program 
will divide tnis number by the total number 
of planet gears in the stage before calculating 
the gear tooth loads . 

2.3. 11 Data Card No. 12 - Shaft Fit Data 


This card is to be included only if the 
change in bearing diametral clearance with 
operating conditions is to be calculated, 
i.e,, if item 3 (ITFIT) on Data Card Ho, 2 
is non- zero . On input data, tight inter- 
ference fits bear a positive sign and 
loose fits a negative sign. 

Item 3 on Card No. 12 is termed the shaft 
effective width. The value specified for 
the effective width may be as large as 
twice the ring width. 
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Use of an effective width is an attempt 
to account for the greater radial rigidity 
of a shaft onto which a ring is pressed. 

This is due to the fact that the shaft 
deflects over a distance that extends beyond 
the ring width. In. the program the 
calculated internal pressure on the ring 
due to its interference fit with the shaft, 
is distributed over the shaft effective 
width. Using double the actual, width as 
the effective width is customary, 

2.3.12 Data Card No. 13 * Shaft Fit Dimensions 
These items are self explanatory. 

2.5.13 Data Card No. 14 ■ Material Properties 

This card defines the modulus of elasticity 
for the shaft, inner ring, rolling element 
and planet gear respectively. This card 
is to be included only if the material 
properties signal on Title Card No, 2 is^l. 

A default value of 204083 N/mm? (29 .6xlO°PSI) 
is used. 

2.3.14 Data Card No. 15 ■ Material Properties 

This card defines the Poisson's ratio for 
the shaft, inner ring, rolling element and 
planet gear, respectively. This card is to 
be included only if the material properties 
signal on Title Card No. 2 is 1. A default 
value of 0.30 is used. 

2.3.15 Data Card No. 16 ■ Material Properties 

This card defines the density for the shaft, 
inner ring, rolling element and planet gear 
respectively. This card is to be included 
only if the material properties signal on 
Title Card No. 2 is 1. A default value of 
7.806 g/cm3 (0,282 lb/in3} is used. 
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Z. 3*16 0ata Card No, 17 - Material Properties 

Tills card defines the coefficient of thermal 
expansion for the shaft, inner ring:, rolling 
element and planet gear respective!/. This 
card is to be included onl/ if the material 
properties signal on Title Card No. 2 is I.. 
A default value of 12.24x10'® 1/®C (6.8x10'® 
1/®F) is used. 


2.3.17 Da ta Card No. 18 - Lubrication and Friction 
DiTt a 

ITEMS 162 

Items 1 and 2 are the amounts by which 
the combined thickness of the lubricant 
film on the rolling track and rolling 
element is increased during the time 
interval between the passage of successive ’ 
rolling elements, frcm whatever replenish" 
ment mechanisms are operative. Item 1 
applies to the outer and Item 2 to the 
inner race"rolling contacts respectively. 

If Item 1 is zero or blank the mode of 
friction is assumed to be dry. 

At the present time the magnitude of the 
inner and outer raceway replenishment 
layers has not been correlated to lubricant 
flow rate, lubricant application methods 
and bearing size and speed factors. The 
user is required to establish proper values 
for the replenishment layer thickness. As 
a rough guide the following suggestions are 
made : 

1) To avoid starvation, the replenish- 
ment layer thicknesses should be 
one to two times the EHD film thick- 
ness which develops in the rolling 
element raceway contacts. 

2) Because of centrifugal force effects, 
intuition suggests that the outer 
raceway replenishment layer should be 
several times thicker than that pre- 
scribed at the inner raceway. 
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ITEM 5 

Item 5, XCAV, describes the percentage of 
the bearing cavity estimated by the user 
to be occupied by the lubricant. 0 ^ XCAV 
< LOO^ This term is used in an equation 
which estimates the amount of frictional 
heat, generated by lubricant churning » 

As with the replenishment layer 'thicknesses, 
the amount of free lubricant should be 
correlated with lubricant flow rate, 
lubricant application methods and bearing 
size and speed factors. At this time such 
correlations do not exist. XCAV values less 
than five percent are recommended. 

ITEMS 4 § 5 


Items 4 and S are the coefficients of 
coulomb friction applicable for the con- 
tact of asperities. . If Items 1 and 2 
are zero, signifying dry operation, then 
items 4 and S serve as the coulomb friction 
coefficients which prevail at respective 
contacts. 

2.3.18 Data Card No. 19 - Lubricant Type 

This card is omitted if Item 1 on card 
18 is zero or blank which implies dry 
friction . 


The relevant lubricant data for four 
specific lubricants (see Table 2) has 
been coded into PLANETSYS. The lubricant 
input information consists only of a 
single number which designates the 
particular lubricant but the relevant 
• information for the lubricant is printed 
in the input data list. 
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TABLE 2 

PROPERTIES OF FOUR LUBRICANTS 
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If Item I, NCODE is 1, 2, 3, or4 PLANETSYS 
uses preprogrsmmed lubricant properties 
' and. no further information is required. 


NCODE 


Lubricant 


I SHELL TURBO 33 

2. MIL-L-7808G 

3" Polyphenyl -Ether 

4 MIL-L- 23699 

NCODE may also be specified as negative 
(-1 to -4) , in which case the traction 
characteristics of the respective 
lubricant NCODE noted above are used 
but the actual properties specified by 
the Items on data cards 20 821 override 
the hard coded data. This option is ^ 

most useful in specifying various mineral 
oils i.e. NCODE = -1. 


If the NASA traction model is being used, two additional lubricant 
data items should be specified on card no. 21: 

a) AKN “ Empirical Lubricant Constant - col. 11-20 

b) FRIG = Lubricant friction coefficient - col. 21-30 


Default values of AKN * 50. and FRIG = .07 will be used if the 
items are left blank. 

For lubricants NCODE ■ 1 through 4 the following values of .4KN 
and FRIG are assigned: 


NCODE 

AKN 

FRIG 

1 

18.2 

,075 

2 

18.2 

.045 

3 

24.9 

.070 

4 

18.2 

.070 
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2.4 Data Set III - Thermal Model Data 

Subsequeat to each, calculation of bearing: generated 
heat rates, either the steady state or transient 
temperature mapping scheme may be executed. What 
follows is a description of input variables, A more 
complete description of the thermal portion of 
”PLANHTSYS’% along with some heat transfer information 
is given in Appendix 

2 .4.1 Data Card No. T1 * System Control Information 

Data Card Tl is a control card. If no 
temperature map is to be calculated, this 
card is to be included as a blank card 
followed by Data Card T2. Data Card Tl 
contains control input for both steady 
state and transient thermal analyses. , 

It is not intended, however, that both 
analyses be executed with the same run. 

ITEM 1 

The highest node number (M) . The highest 
node number must not exceed one hundred, (100). 

ITEM 2 

Number of unkno^m temperature nodes (N) . This 
number should equal the total number of unknown 
temperatures. It is required that all nodes 
with unkno'im temperatures be assigned the 
lowest consecutive numbers starting from one 
(1). The program assumes that all node numbers 
greater than N (from N+1 to M) represent known 
boundary temperatures. 

ITEM 3 

Common Initial Temperature (TEMP)°C: The 

temperature solution iteration scheme requires a 
starting point, i.e., guesses of the equilibrium 
temperatures. Data Card T3 allows the user to 
input guesses of individual node temperatures. 
When a node is not given a specific initial 
temperature, the temperature specified as Item 
3 of Data , Card Tl is assigned, 

24 
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ITEM 4 

Puncli Flag CIPUNCH): . If the Punch Flag is not 
xero (0) or blank, the system equilibrium 
temperatures along with the respective node 
numbers will be punched according to the format 
of Data Card T3. This option is useful if the 
user makes a steady state run with lubricatiori 
and then wishes to use the resultant temperatures 
for initialization of a transient dry friction 
run* The latter is performed to assess the con- 
sequence of lubricant flow termination* 

ITEM 5 

’’Output Flag" (lUB) . If the "Output Flag" is 
not zero the bearing program output and a 
temperature map will be printed after each 
call to the planet bearing solution scheme. 

This printout will allow the user to 
observe the flow of the solution and to note 
the interactive effects of system temperatures 
and bearing heat generation rates. Since the 
temperature solution is not mathematically 
coupled to the bearing solution the possibility 
exists that the solution may diverge or oscillate. 
In such a case, study of the intermediate output 
produced by the "Output Flag" option may provide 
the user with better initial temperature guesses 
that will effect a steady state solution. 

ITEM 6 

"Maximum Number of Calls to the planet bearing 
Program" (ITl) . ITl is the limit on the 
number of Thermal-planet-bearing iterations. 

The user must input a non-zero integer such as S 
or 10 for the Program to iterate to an equilibrium 
condition. If ITl is left blank or set to zero 
(0) or 1, planet bearing performance will be 
based on the initially guessed temperatures of 
the system. The temperatures printed out will 
be based on the bearing generated heats. It is 
unlikely that an acceptable equilibrium condition 
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will be achieved. However, the temperatures 
which result may provide better initial guesses 
for a subsequent run than those specified by 
the user, 

ITl also serves as a limit on the transient 
temperature solution scheme, by limiting the 
number of times the planet bearing solution 
scheme is called. Each call to the 
planet bearing scheme will input a new set 
of bearing heats to the transient temperature 
scheme until a steady state condition is 
approached or until the allowed transient 
solution time limit is reached* 

ITEM 7 

'^Absolute Accuracy of Temperatures for the 
External Solution" (EPl) , In the steady 
state thermal solution scheme, each calculation 
of system temperatures occurs after a call 

to the Diane t- bearing routine which produces 
bearing generated heats. After the system 
temperatures have been calculated for each 
iteration, using the temperature mapping 
heat dissipation scheme, each node temperature 
is checked against the nodal temperature at 
the previous iteration. 

If {^n,i “ ^n-l,i} _< EPl for all nodes i, then 
equilibrium has been achieved and the iteration 
process stops. Where ^n, i =! temperature of 
ith node at nth thermal iteration; ^n-i, i = 
temperature of ith node at n-lth thermal iteration. 

ITEM 8 

"Iteration limit for the temperature mapping 
heat dissipation Solution" (IT2). After 
each call to the planet bearing program, the 
node to node heat transfer scheme is solved 
to determine the steady state equilibrium 
temperatures throughout the system, based on 
the calculated set of bearing heat generation 
rates . If IT2 is left blank or set to zero, 
the number of iterations permitted to solve 
the heat transfer scheme is limited to twenty, 

( 20 ). 
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ITEM 9 

•’Accuracy for Internal Thermal Solution" (EP2) . 

The use of EP2 is explained in Appendix B (page 
Ell), If EP2 is left blank or set to zero (0), 
a default value of 0.001 is used. 

ITEM 10 

"Starting Time’* (START) is a time, T_, at 
which the transient solution begins; usually 
set to zero (0). 

ITEM 11 

"Stopping Time" (STOP) is the time in seconds 
at which the transient solution terminates, 

T,. The transient solution will generate a 
history of the system performance which will 
encompass a total elapsed time of 

(Tj - Tg) seconds. 

ITEM 12 

"Calculation Time Step" (STEPIN). The transient 
internal solution scheme solves the system of 
equations: 


^k+1 - + % At/pCpV 

ht - STEPIN 

The user may specify STEPIN. If left blank or 
set to zero (0), the Program calculates an 
appropriate value for STEPIN using the procedure 
described in Appendix E (pp. Ell -ElS). 

ITEM 15 

"Time Interval between Printed Temperature 
Maps*' (TTlMB) seconds. The user must 
specify the length of time which will 
elapse between each printing of the temperature 
map. The interval Av’ill always be at least 
as large as the "calculation timestep" 

(STEPIN) . 

2i 
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ITEM 14 

"Time Interval Between Calls of the Planet 
Bearing Portion of the Program" (BTIME) . 

BTIME will always have a value larger than or 
equal to (STEPIN) even if the user inadvertently 
inputs a shorter intervals Computational time 
savings result if BTIME is greater than STEPIN; 
however), accuracy might be lost. 

2.4.2 Data Card T2 - Bearing Temperatures 

Card Type T2 is required, one card for each 
planetary stage, if no thermal analysis is 
being performed. The « temperature data is 
used within the planet-bearing analysis to 
fix temperature dependent properties of the . 
lubricant. 

For this purpose the inner race, outer race 
and lubricant bulk cavity temperatures are 
used in the analysis which calculates the 
change in bearing diametral clearance from 
"off the shelf" to operating conditions. 
Temperatures must be given in °C, 

2.4.3 Data Card T5 - Nodal Temperatures 

In the steady state analysis this card is 
used to input initial guesses of individual 
nodal temperatures for unknown nodes as well 
as the constant temperatures for know nodes, 
such as ambient air and/or an oil sump. 

In the transient analysis, Data Card T3 
is used to input the nodal temperatures 
of all nodes at (START) » Ts, i.e., at 
the initiation of the transient solution. 

2 • • 4 Data Card T4 - Bearing Node Numbers 

With this card, node numbers are assigned 
to the components of each bearing, one card 
per stage. With this information the proper 
system temperatures are carried into each 
respective bearing component to be used in 
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the bearing analysis. The inner race and ihhejc ring 
node numbers may or may not be the same at the user's 
discretion. Similarly the outer race and outer ring 
node numbers may or may not be the same, 

2.4.5 Data Card T5 - Heat Generating Nodes 

Card type T5 is required, one card for each planetary 
stage, if a thermal analysis is to be performed. The 
planet bearing system analysis accounts for frictional 
heat generated at five locations in the bearing, i.e. 
at the inner race, the outer race, between the cage 
rail and ring land, at the cage rolling element con« 
tacts and in the bulk lubricant due to drag. This 
card allows the heat generated to be distributed 
equally to two nodes. For instance, the heat generated 
at the inner race-rolling element contact should be 
distributed half to the rolling element and half to the 
inner race. The heat developed between the cage and 
the inner ring land may be distributed half to the in- 
ner ring and half to the cage, if a cage node has been 
defined, otherwise, half to the bulk lubricant. 

2.4.6 Data Card T6 - Nodes with Constant Power Source 


This card specifies the node numbers and the heat 
generation rate for those nodes where heat is generated 
at a constant rate such as at rubbing seals . 


The heat generated by the planet-sun and planet-ring 
gear contacts is not currently calculated by PLANETSYS . 
In order to do a thermal analysis of a planetary bear- 
ing assembly, the heat generation rates for these con- 
tacts must be calculated by the user and input using 
T6 data card (S) . The heat distributed to the plane- 
tary gear node should correspond to a single planet 
contact. The heat distributed to the non-planetary 
nodes (sun gear, ring gear) should correspond to the 
total for all planet contacts. 


2.4.7 Data Card T7 - Heat Transfer Coefficients 

This card type is used to input the numerical values 
of the various heat transfer coefficients which appear 
in the equations for heat transfer by conductivity, 
free convection, forced convection, radiation and 
fluid flow. Up to ten coefficients of each type may 
be used. Separate values of each type of coefficient 
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are i^ssigned an index number via card T7 and in 
describing heat flow paths (Data Card T8 below) it is 
necessary only to list the index number by which heat 
transfers between node pairs. 

Indices 1-10 are reserved for conduction, 11-20 for 
free convection, 21-30 for forced convection, 31-40 
for radiation and 41-50 for fluid flow (product of 
specific heat, density and volume flow rate) . 

As an example, for heat transfer by conduction within 
steel having a thermal conductivity of 53.7 watts/m® c, 
one could prepare a card type T7 with the digit 1 
punched in column 10 and the value 53,7 punched in the 
field corresponding to card columns 11-20. If a con- 
duction coefficient of 46.7 were applicable for certain 
other nodes in the system one could punch an additional 
card assigning index nq. 2 to the value 46,7 by punch-' 
ing a "2” in card column 10 and 46.7 anywhere within 
card columns 11-20 . 

Rather than inputting constant forced convection 
coefficients, optionally, these coefficients can be 
calculated by the program in one of three ways, if 
the calculated option is exercised, a pair of cards is 
used in place of a single card containing a fixed value 
of a. The contents of the pair of cards depends upon 
which of the three optional methods are used. 

Option l)a is independent of temperature, calculated 
as a function of the Nusselet humber which 
in turn is a function of the Reynolds num- 
ber Re, the Prandtl number Pr as 
follows, (cf . 17) . 

« = Nu X oil/L 
Nu= aRe^Prc 

Where A is the lubricant conductivity, L is 

a characteristic length (with a unit of meters) 
and a, b, and c are constants. 
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Option 2) a is a function only of fluid 

dynamic viscosity, and viscosity 
is tampetature dependent. 

a » cr|.<^ 

where c and d are constants 

Option 3) a is a function of the Nusselt, 
Reynolds and Prandtl numbers and 
viscosity is temperatur® dependent . 


2.. 4. 8 Data Card T8 - Heat Flow Paths 

This card defines the heat flow paths between 
pairs of nodes. Every node must be connected 
to at least one other node, i.e., two or more 
independent node systems may not be solved 
with a single Program execution. 

The calculation of heat transfer areas is 
based on lengths, L, and input using Card 
Type 8. Additionally, the type of surface 
for which the area is being calculated is 
indicated by the sign assigned to the heat 
transfer coefficient index. If the surface 
is cylindrical or circular the index should 
be positive, if the surface is rectangular the 
index should be input as a negative integer. 
Note that due to the axial symmetry of the 
planetary system most nodes will represent 
circular rings and thus the heat transfer 
indices are predominantly positive. For 
conduction the length L- gives the separation 
distance between the two nodes. 

In the case of radiation between concentric 
axially symmetric bodies, L- is the radius 
of the larger body. For radiation between 
two parallel flat surfaces or for conduction 
between nodes, L- is the distance between them. 
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Fluid flow heat transfer accounts for the energy which 
the fluid transports across a node boundary. Along a 
fluid node at which convection is taking place, the 
temperature varies. The nodal temperature which is 
output is the average of the fluid temperature at the 
output and input boundaries. If the emerging tempera- 
ture of the fluid is of interest, it is necessary to 
have a fluid node at the fluid outlet . At this auxil- 
iary node only fluid flow heat transfer occurs and the 
fluid temperature would be constant throughout the 
node. Thus the true fluid outlet temperature will be 
obtained. 

Conduction of heat through a bearing is controlled by 
index 51, The actual heat transfer coefficients which 
contain a conductivity, area and a path length term 
are calculated in the bearing portion of the program. 
The term is based upon an average outer race and inner , 
race rolling element contact. 


Special guidelines must be followed in specifying heat 
•transfer linkages between planetary nodes and non- 
planetary nodes (e.g., conduction between planet gear 
node and ring gear node) . For conduction, convection, 
or radiation heat transfer, the planetary node must be 
defined as node I (columns 11-15) . The number of 
planets must be input in columns 51-60 . The heat trans- 
fer area defined by L.^ (columns 21-30) and L 2 (columns 
31-40) must correspond to a single planet. 

For fluid flow heat transfer between a planetary node 
and a non-planetary node, the data on card T8 depends 
on the direction of flow. The program assumes that the 
direction of fluid flow is from node I (columns 11-15) to 
node J (columns 16-20) . When node I is the planet 
node, the first fluid flow index (columns 1-10) must 
correspond to the flow from an individual planet. The 
second flow index (columns 21-30) must correspond to 
the total flow from all planets. The number of planets 
is specified as a positive number in columns 51-60 . 

When node J is the planet node, the first fluid flow 
index must correspond to the total flow to all planets. 
The second index must correspond to the flow to an in- 
dividual planet. The number of planets is specified as 
a negative number in columns 51-60 . 
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9 Data Card T9 * Node Heat Capacity 

This card Inputs data required to calculate the 
heat capacity of each node in the system. This 
card type is required only for a transient 
analysis. 

3 . COMPUTER PROGRANt OUTPUT 
3.i Introduction 


The Program Output is intended to provide the user 
with a complete picture of the planet bearing 
system performance. 

This portion of the Manual is not meant to be a 
list and explanation of all program output. Most 
items are self explanatory and have been omitted 
from further discussion. 


Two sample program outputs are included in Appendices 
B and Cji reflecting use of the NASA and SKF traction 
models, respectively. In both cases, the output 
represents the solution for a single stage planetary 
gear system transmitting 202.2 HP. 


The first five pages of output essentially consist 
of a summary of the input data categorized into 
bearing, cage, steel, lubricant and loading data. 

The remaining six pages contain the calculated program 
output. 

3 . 2 Bearing Output 

3.2.1 Gear Tooth and Inertial Loads 

Figure 4 shows a single gear tooth and 
the applied loads. The gear tooth torque 
arises as a result of the tangential tooth 
load with a lever arm equal to the difference 
in radii between the planet pitch circle 
and planet ring neutral axis. 
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FIGURE 4 IDEALIZED GEAR TOOTH LOADS 
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3.2,2. Fatigue Lives 

Tixe fatigue life of the outer and 
inner‘s raceways as well as the bearing are 
calculated based upon the 4th power relationship. 

The; bearing life represents the statistical 
combination of the two raceway lives and in 
the case of a double row bearing, the combination 
o£ the two single row lives which are assumed 
to be equal. These lives reflect the combined 
effects of the lubricant film thickness and 
material life factors. 


3.2.5 h/sigma 


The ratio h/o also referred to as A, , is 
printed for the most heavily loaded rolling 
element. The variable h, represents the 
EHD plateau film thickness with thermal 
and starvation effects considered. The 
variable cr, represents the composite root 
mean square surface roughness of the rolling 
element and the relevant raceway. 

3.2.4 Life Multpliers 

Lubrication - This life multiplier is a 
function ot h/a at each concentrated 
contact. Its value ranges from 0.21 for 
A ^ 0.6 to 3.0 for A ^ This subject 

is covered in more detail in Section 2.3.3, 
Item 5. 

Material - This output simply reflects the 
input ue . 

Again, it is covered in Section 2.3.3. 
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temperatures Relevant to Bearing Performance 

These temperatures fully describe the 
temperature conditions which affect the 
performance of a given bearing. If one of 
the temperature mapping options is used, 
the temperatures printed reflect the 
resiLtts of the particular option. If 
neither temperature option is used, the list 
is simply a repeat of the input data. Note 
that there are separate temperatures for 
outer and inner raceway and ring temperatures. 
The raceway temperature is used to determine 
lubricant properties. The ring temperatures 
are used in the bearing dimension change 
analysis. The raceway and ring temperatures 
may be the same value. 

3.2.6 Frictional Heat Generation Rate and 
flearing Friction Torque 

Frictional Heat Generation Rate 

The various sources of frictional heat 
generated within the bearing are listed. 

The values printed for "OUTER RACE, INNER 
RACE, R.E. DRAG, R.E. CAGE and CAGE- LAND" 
represent the sum of the generated heats 
for all rolling elements and cage, Add- 
itionally, the heats printed for the 
outer and inner raceways plus the rolling 
element cage, reflect the friction 
developed outside the concentrated contacts, 
i.e. the HD friction^ as well as the EHD 
friction developed within the concentrated 
contacts. "R.E, DRAG" should be interpreted 
as the heat resulting from lubricant churning 
as the rolling elements plow through the air- 
oil mixture. Cage-land heat generation rate is 
an estimate since unavoidable variations in 
cage-land diametral clearance have a marked 
effect on cage-land temperature, making precise 
calculations impractical . (See ref. 15) . 
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Torque 

The torque value i 5 calculated as a function 
of the total generated heat and the suin of the 
outer and inner ring rotational speeds. The 
intent is to present a realistic value of the 
torque required to drive the bearing. 

S *2 . 7 EHD Film and Heat Transfer Data 

EHD Film Thickness 

These values refer to the calculated EHD plateau 
film thickness at both contacts of the most 
heavily loaded rolling element and include 
the effects of the thermal and starvation 
reduction factors. 

Starvation Reduction Factor 

These factors give for the inner and outer 
ring contacts, the reduction in EHD film 
thickness due to lubricant film starvation 
according to the methods of Chiu, (9). 

These factors pertain to the EHD film 
thickness for both the inner and outer race 
contacts of the most heavily loaded 
rolling elements, but are applied to the 
respective inner and outer race film thickness 
for each rolling element in the bearing, 

Thermal Reduction Factor 

These factors are calculated according to the 
methods of Cheng, C8) and pertain to the EHD 

film thickness for both the inner and outer 
race contacts of the most heavily loaded 
rolling elements, but are applied to the 
respective inner and outer race film thickness 
for each rolling element in the bearing. 
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Meniscus Distance 


This parameter is defined as the distance 

along- the direction of rolling from the 

contact center to the oil meniscus in the 

contact inlet. It is calculated according to 

the methods of Chiu C9l for the inner and outer 

race contacts of the most heavily loaded rolling element 

Raceway- Rolling Element Conductivity 

These data reflect the amount of heat 
transfer between rolling element and 
raceway per unit degree difference between 
the two components. These data reflect the 
average of all outer and inner contacts 
respectively. ‘ 

3.2.8 git and Dimensional Change Data 

Fit Pressures 

These data refer to the pressures built 
up as a consequence of interference fits 
between shaft and inner ring and housing and 
outer ring. Pressures are presented both 
for the standard cold-static condition (16° C) 
and the operating condition. 

Clearances 

"Original" refers to cold unmounted clearance which 
is specified at input if the diametral clearance 
change analysis is executed. ■ "Change" refers to 
the change in diametral clearance at operating 
conditions relative to the cold unmounted 
condition. A minus sign indicates a decrease 
in clearance . "Operating" refers to the clearance 
at operating conditions. 
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» 1*1.9 Lubricant; Temperatures and Physical Properties 

The- lubricant properties, particularly the 
^ dynamic viscosity and to a lesser degree the 

pressure viscosity coefficient are heavily 
* temperature dependent. These factors enter the 

EHD film thickness calculation and the HD and 
-EHD friction models., The lubricant is assumed 
ii< to be at the same temperature as the relevant 

raceway. As noted elsewhere, these temperatures 
f* may be either input directly or calculated by 

y the Program. 


The physical properties printed are self 
explanatory. The units are enumerated. 

3 • 2 . 1 0 Cage Data 

Cage -Land Interface 


[ 

[ 

[ 

C 

I 

15 


^ i 


The cage. data indicates the performance parameters 
at the interface between the cage rail and the 
ring land on which the cage is guided. The 
torque and heat rate require no explanation. The 
eccentricity ratio defines the degree to which 
the cage approaches the ring on which it is 
guided at the point of nearest approach. The 
radial displacement of the cage relative to the 
bearing axis is divided by one half the cage-land 
diametrical clearance. An eccentricity ratio of 
one indicates cage-land contact. A ratio of 
zero indicates that the cage rotation is con- 
centric with the bearing axis. The cage -land 
separating force results from the rolling element 
cage loads driving the cage to an eccentric 
position and causing pressure to develop between 
the cage rail and ring land. Only the cage -land 
and rolling element pocket forces are Gonsidered 
in determining the cage eccentricity. The 
centrifugal force which results from the 
eccentricity, although available, is not 
considered in the analysis. 
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5. 3 Rolling Element Output 

S,3«L Rolling Element: Orientation 

The first rolling element may be located 
in either *o£ the following ways: 

■ li. When the planet bearing center is 
"■ stationary as in cases 1 and 2 of 
Table 1, the first rolling element 
is encountered at 9 o'clock, gear 
tooth loads being applied at 12 and 
6 0 ' clock. 

2. When the planet bearing center 
orbits as in cases 3 thru 6 of 
Table 1, the first rolling element 
is located at 12 o’clock plus half 
the rolling element spacing. Again, 
gear tooth loads are applied at 12 
and 6 o’clock. 

3.3.2 Rolling Element Raceway Loading 
Normal Forces 

The gear tooth loads, the inertial load of 
the planet gear, and the rolling-element 
contact loads acting on the planet bearing 
outer ring form a self-equilibrating loading 
system. 

The rolling-element contact loads depend on 
the bearing outer ring deflections, the 
rolling-element inertial loads, and the radial 
displacement of the bearing inner ring, which 
is considered as an elastic solid. 

The cage force is always directed tangent to 

the bearing pitch circle. If the rolling 
element orbital speed is positive, a positive 
cage force indicates that the cage is pushing 
the rolling element, tending to accelerate 
it. Cage force is a function of the tangential 
component of the rolling element inertia less 
the raceway friction forces. 
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Hertz Stress 

The stress printed represents the maximum 
normal stress at the most heavily loaded 
slice of tlie roller raceway contact. 

Outer Ring Deflections 

For a flexible outer ring, the ‘radial 
deflections are given at the various rolling 
element locations. The sign convention of 
radial displacement is shown in Figure S. 

5.3.3 Deflection Plotting 

A plot showing the outer ring deflections 
and rolling element - outer race contact 
loads versus azimuth angle is provided at 
the User's Option. (See sect, 2.2,2, Item 6) 

3.4 Thermal Output 

As is the case for bearing output, all of the 
input data is printed. The calculated output 
data is presented in the form of a temperature 
map in which a node number and the respective 
nodal temperature appear. The appearance of the 
steady state and transient temperature maps are 
identical.. The transient temperature map also 
includes the time (T) at which the temperature 
calculations were made. 

4 , GUIDES AND LIMITATIONS TO PROGRAM USE 

''PLANETSYS" is intended to be as general as possible 
with the following limits on system size: 

1. Number of rolling elements per bearing 
row must not exceed 40, 

2. Gear tooth loads must act at the two 
ends of the diameter on a line passing 
through the sun and planet gear centers 
and in the radial plane containing the 
bearing center. 
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3. Number of stages in the transmission 
system is limited to 3, 

4r Number of temperature nodes used to 
describe the system - 100 maximum, 

> Some general guides for efficient use of the Program 
are listed below: 

1. Attempt to input bearing operating 
diametral clearance rather than cal c ul a t e 
it. Of, execute the diametral clearance 
change analysis once for a group of similar 
runs and use the output from the first run 
as input to the subsequent runs omitting the 
clearance change analysis, 

* 

2. Attempt to input accurate operating 
temperatures rather than calculate them. , 

3 . The more non-i inear the problem the 
more computer time . In the the rma I' 
solutions, if possible, eliminate 
nonlinearities by omitting radiation 
terms and by using constant rather than 
temperature dependent free and forced convection 
coefficients. 

4. In the transient thermal solution, space 
the calls to the planet-bearing solution 
TFTimej to as large an interval as prudently 
possible . Be careful, however, too long 

an interval will produce large errors in 
heat rate predictions. 

5. In the steady state thermal analysis at temp t . 
to estimate nodal temperatures on a node by 
node basis. Nodes which are heat sources 
should have higher temperatures than the 
surrounding nodes. 
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It is also suggested, that a constant user of the ptogram 
should stud/ the hierarchical Program flow chart Appendix 0, 
alony with the Program listing to gain an appreciation 
of the program complexity and the flow of the problem 
solution. The Program is comprised of many small functional 
subroutines. Knowledge o£ these small elements may allow 
the user to more easily piece together the philosophy of 
the total problem solution. 
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APPENDIX A 

»‘PLANEtSYS»* INPUT DATA FORMS 
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Ho. of Hedtiction 
Stages - 3 Max. 


Debug yiag 
IF 0 or Blank 

No Debug Output 

IP - 1 Bain Loop 
Debug Only 

IF * 2 Full Debug 
Output 


Pit Calc Flag 
0 or Blank 

No Fit calcs — 
and no 
iteration* 

Any positive No- 
Fit calcs perf, 
witb a xaax of 5 
' iteration* 

No. 

Fit calc* are 
perforaed. No. 
of iteration* 
equal* ABS • 
value of given 
nuiRber 


Main Loop Iter- 
ation Flag 

Any positive no. 
Fre-Set Max 
(20) i* used 

No- of iter- 
ations equals 
ABS. vetlue of 
given number 

0 or Blank 
One iteration 


Fit Loop Accuracy* 

0 or Blank 
Pre-Set 
accuracy of 
.0001 is used 


Plot Flag 

Blank or ^False^ 

Eo Plots 

■True* 

Cutout will 
contain plots 
of Outer Ping 
deflection v*. 
azinuth angle in 
cartesian 
coordinates 


Dimensions Flag 

0 or Blank 
Input variables 
are given in 
English 0nits 

1 Input variables 
given in Si 
Units 


Material Prc^rfcie 
Signal 

^-0«er will 
Specify «at. 
properties 

0 or Blank 

Pre-set material 
properties are 
used omit data 
cards 14 through 



♦ Fit Loop Accuracy defined as change in dla^trical ^^learance 
divided by Rolling Element Diameter. 
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DATA CARD 18 - Cage Xnfonaation 


BD(44) 


BD(45) 


n 1 0 Jo 


Cage Type 


F I 0 40 


Rail-Land Diaiseter 


BD(46) 


FI 0 . 0 




BD{47) 




Ball width 


32^3 


F I 0 JO 




3eb3 


BD(48) 




1 0 


Ball Land Oia. 
Clearance 
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{clearance 
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DATA CAiU> 110 
BD(32) 


BD(33) 


BD(38) 


BD{39) 


BD(40) 
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BD(42) 
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No* of Planetary 
Bearings in Stage 



ee Table 1 
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Planet Geeur 
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Width of Planet 
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Planet Outer Ring. 
Weight 
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DATA CARD #21 - Cube Properties 

BDC98) 


0*1 1 this card If outer raceway replenlshaent. lay< 
.data card #18, is 0 or blank. Also omit the card 
data card 119, is 1-**- Bead this data If HCODE < 


, \.h Ui. UIt lel9 |.ol» |>2|u|i4hsl.6|i7|ie|iaN2»N23|^^NN^^^^ 


FI 0.0 


FI 0.0 


THERMAL 

CONDUeTIVITV 


BTU/FT - r®-HR 
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AKN 

EIID niGH 

CONTACT-STRESS 

FACTOR 


F 1 Oi 


FRIC 

ALLEN 

FRICTION 

COEFFICIENT 


INPUT THESE VALUES VIIIEN 
EXECUTING NASA VERSION 
(SEE PAGE 23 AND APPENDIX F) 
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BATA CAHO T4 — BEARING NODE NUMBERS Include one T4 C 2 urd for each planetaxy sta^e*^ 
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^2 
TI o 
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O tj 
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jiata CAivn t*i 


lit taifiutATion koih-5 


ttsr inittt* out* rS car4l «ac;ii {ilaiteicsiry sta^e^ 
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I 
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oo 
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> 
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OO 
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bATA CA«n T6 


thijfm;)! with a i>ower aoutHHi 

UtlB card.% nena^d . followc-d by a blank card 
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^SPER COEFFICIENTS Input one or two cards per coefficient. Use as inany T7 cards as required followed by a blank car . 


iiiiiiiiMHiiiiiiiBMiiiMmMMiiiiiiiiiiiiiiiiHiHmiii 


INDEX FOR LATER 
IDENTIFICATION OF 
HEAT TRANSFER 
COEFF. 

1-10 ( Conduction) 
11-20 (Free Conv. ) 
21-30 (Foroed Oonvt) 
31-40 (Radiation) 
41-50 (Mass 

Transport) 


MAGNITUDE OF HEAT EXPONERT USED WITH 
TRANSFER COEF- TEMP. DIFFERENCE 
FICIENT 


l(w/m-Oc) 

ay(w/ra2-Oca) 

a^w/m^-Oc) 

’Lo^Cw/Oe) 


a (default 


(Q = AA AT/1) 

(Q = avA AT*) 

(Q^- awA AT)^ ^ 

p = BensTty (KgJm^), Cp = Specific Heat (W-SEC/KG-Oc) , v = Flow Rate Cm3/SEC) 


•The forced convection coefficient can be calculated Internally by the program using one of three options. For options 1 and 3^ a^, 
is calculated by ay ^ X Nu/L where Nu = KReApr® (Re = Reynold's No. = ULp/n, Pr = Prandtl No. = nCp/X). For option 2, ay = ctj^- 
The viscosity is a function of temperature for options 2 and 3* To activate this feature, input data for card T7 as shown below 
and follcw immediately with card T7A. 


HEAT TRANS. INDEX 
21-30 
21-30 
21-30 



U (M/SEC) 
U (M/SEC) 


X (W/m-*C)^ 
X (W/m-*C) 


^2 

-n 2 

*tl Q 
OS 
o > 

50 p 

o -o 


r " m 


(OPTIOK 1) 
(OPTION 2) 
(OPTION 3) 


DATA CARD T7A Input only xf it 13 desired to calculate forced convection coefficient internally by 
the program. 


iigolBasliii§^iaiiiiii^iiisiii3iiiiiEBEsaBS3SS3asaxsiasss3sssss3£S£iai§2aasaaas 

niiiiiiiiiiimiiiiiiiiiiiiimiiiniiiiiiimiHiiiiiiiHiiiiiiiiHimr — 


yr(M-’Secy>l^) 

p(KG/H^) 

Cp(W-sec/KG-®c) 

— 

p (KG/M^) 

Cp (W-scc/KG-®c) 


fmmm jmmh 



n § T^ (ll-sec/M"^) 

n 8 T_ (N-sec/M^) 
Ij 





ro 


Ccy 



(OPTION 1) 
(OPTION 2) 

(OPTION 3) 




















data CAAD T8 

Heat Flow Paths 

USE As Many Cards As Needed, Max, 500, Followed By A Blank Card 


[IIS3fl3SEIQSBBIS9&BQI3SSSSSSSSSSSSS9SSSSSISSSSS3ESS9SSSS9SS 


INDEX = 51 


NODE I 

NODE J 

NODE I 

NODE J 

NODE I 

NODE J 

NODE X 

node j 

NODE I 

NODE J 

NODE I 

NODE J 

^ODE I 

. 

NODE J 


Index of fluid 
Flow at NODE J, 
41 < INDEX <50 


1 ^ NO. < 3 


I»2 im) 


1^2 (bob) 


Raceway Flag 

1. Inner Race 
Contact 

2. Outer Race 
Contact 


No- of planets in the stage. Input only if .heat 
transfer involves a planetary subsystes node and a 
non-planetary node. See section 2.4.8 for detailed * 
instructions* . ^ 


Conduction between I and J area = 2^ 

If index < 0 area — Distamce I-J = 


Natural convection between I and J. 

Area = 2ir I* 2 l‘ 2 * index < 0 area »= ^x^2* 


Forced convection between I and J. ^ 
Area as above. If n(t) , t is t . 


Radiation between J and d. Area as above. 
For description of see User's Manual. 


I to node d. First index is index 
flow at node I. Second index corresponds to 
Jfluid flow from I to J. 


O -O 
c > 
> o 
cm 


o 

o 
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APPENDIX B 
S.AMPLE OUTPUT 
"NASA” VERSION 



SKF TECHNOLOGY SERV'CES 

SKFlNDUSTOlES.iNC 


B-1 



P L A H E T S Y S •••• technology DIVISION S K F INDUSTRIES INC. •••• P W * N f T ST S 


planctsys executed using n as a traction kooea 


THIS planetary SYSTEH CONTAINS I STAGE. 


STAGE NO, (l> CONTAINS 3 IDENTICAL PLANETS, EACH SUPPORTED BY A SPHERICAL ROLLER BCABIH® HAVIJli 2 ROlfS OF ROLLERS. 
THE MAXIHUI HO. OF MAIN LOOP ITTERATIONS ALLOWED IS 20 AND THE RELATIVE ACCURACY REQUIRED IS ,10000-07 


w 

i 

to 




AT81D044 




I’ I. 4 M r T S T S TECHNOlOGT division S K F IHOUSTRILS JHC 


P L A N E T S V S 




I LA.'it rSYS EXECUTED USiriG H a' S A TRACTION MODEL 


I 


Cd 

I 

ot 


! 

i 

1 

fe 


N U T E S nDIKCNSiONAL UNIT^ ANGLES .......................... .DECREES 

density*.. ............. POUNOS/CUBIC INCH 

elastic MODULI ... ...... ..POUNOS/SQ. INCH 

FORCE...... ....POUNDS 

KINEMATIC Vise .......... .inches SO. /SEC. 

LENGTH.... ......INCHES 

HASS........ .POUNDS 

PRESSURE POUNOS/SQ. INCH 

SPEEDS. ....R.P.M. 

SURFACE ROUGHNESS. ...........MICROINChES 

TEMPERATURE. ...... ....DEGREES FAHRENHEIT 


THERM. CONO........B.T.U./HR-rOOT-OEG F. 

THE term bearing number WHICH APPEARS BELOW IS SYNONOMOUS WITH STAGE NUMBER • 

3) planetary inversion refers to The specific kinematic condition as oescribco in the users hanual. 


> 

H 

00 


o 


. RAmiliMiiiN Vnii.iiiii 
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B-4 


* 1 A n r T S Y S TECH'JOLOCY uivision s k f 

PLAhCTSYS EXeCUTEU USIHtf N ASA IHACTION lODEL 


nmaea 

I 

pitch dm 

Or RCARirJi: 

.2A»01 

• * » 

PLANET 

MASS 

3.4A0 

^planetar 

PITCH DIA 
OF GEAR 
.395*01 

BrARI^/f; 

«u«be*i 

1 

Nu^iaci 
OF PLAHETS 
3 

POWER THSU 
STAGE 
.202.200 

plan, invers. 
ISEE NOTE 3i 
3 


INDUSTRIES 

p L A 

>• E T S T S 


Y OAT A 
WIDTH OF 
PLANET 
1.550 

0,1 A O F 

PL. NTRL. AXIS 
3.256 

PLANer HONCflT 
OF IHEgirA 
fOU 

planet TOOTH 
PRESS, angle 

2A.6 

SPEED OF 
SUN 

1659.Q 

SPEED OF 

carrier 

355.0 

SPEED OF 
PING 

• *0 

PITCH OlA. 
OF SUN 

3.0 
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. P L A « C T S Y S •••• TECIIMOLOOY OIVISION S K F IMDUSTRIES IltIC. •••• PUAI(E1STS 
PUAHrT'iYS EXECUTCO USIHG A 3 * TRACT !OM HOOEL 


llEAklM'5 

NUnijE.); 


* < * • • 8 C A R I N S 
BEAR I NT, type 

SPHERICAL roller BEARING 


DATA 
NO. OF 
ROUS 
2 


NO. OF R0LLIN6 CONfACT 
CLTS. PER ROU AHCLE 

12 IS.fl 


OIAHETRAL 

C|.EARANCE 

.00230 


bEiiU INb 
mHtiZH 
I 


cage: ttpe 

ROLLING LLCHEMT RIOING 


: • C A C E D A I A « « • 

CAGE POCKET RAIL LAND 

CLEARANCE UIOTH 

•0100 .0000 


RAIL LAND RAIL L#NO 

oianeter clearance 

•000 •OOi 


CAGE 

HASS 
• OA 


Cd 

I 

^ aCARrriG 
MUrirtEK 

1 


...steel data. • • 
steel type. life steel type* 

INKER RIN3 FACTOR OUTER RING 

H-nn FOR-7 1.000 CARBURIZED 9310 


LIFE 

FACTOR 

1.000 


r 

i 


> 
>-3 
oo 
M 

a 

o 

s ■ . 
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.> P L A M C T $ Y S technology DIVISION S K F INDUSTRIES INC* »••• P L A N ( ^ S Y S 

PLAMETSYS EXECUTED USING NASA TRACTION BODEL 


ROLLING ELEBENT DATA 


OEARING NO.tUi SPMCR ICAL .ROLLER HE AR INO 

ROLLER OIAIETER ROLLER TOTAL LENGTH ROLLE.RPROFILE RADIUS NO» OF AX|AE SLICES 


.5118 

.5520 

I.A5A0 

II 

OUTER RACE WAY 
EFFECTIVE LENGTH 

OSCULATION 

INNER RAtCIIAT 
EFFECTIVE LENGTH 

OSCULATION 

.A520 

• 9690 

.9520 

*9690 ■ 





ORIGINAL PAGE IS 
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* P t A fl e T S Y S •••• TECHNOLOGY DIVISION S K F IWOUSTRIES INC. •••• PLANC1SYS 
PLA'JL'TSYS EiCECUTEO USING N A S A TRACTION MODEL 


. . . . . S Ij R F A C E D A T A • • • • • 

UEARlMa CLA ROUOHNESS RKS ASPERITY SLOPE 

NUHDER outer inner roll. elm. outer inner ROLL. CLN. 

L 8.00 8.00 6.00 l.OOO 1.000 1.000 


• * • . . L U 0 a I C A NT 0 A T A • • • • • 

UEARINJ DESIGNATION KINEMATIC VISCOSITY DENSTTY AT THERMAL EXPAN. THERMAL 

NUnUER (37. 7R C) (98.89 C> (15.56 C) COEFFICIENT CONDUCTIVITY 

I HIL-L-23699 .04 .01 .04 . 41-03 .088 


. 4*. .lubrication AND FRICTION DATA* 

® BEARING PCRCE'IT LIIOE FILM REPLENISHMENT ASPERITY 

^ JJ'JHUER IN CAVITY LAYER THICKNESS FRICTION 

(ROLL.ELM. ♦ RACEUAY I COEFFICIENT 

OUTER INNER 

1 1.00 .4000-04 .2000'‘04 .10 

GIVEN TEHPERATURES ( OCG C > 

*m; Shaft i. ring i. race roll el. o.race planet cage bulk 

1 121.00 121.00 121.00 121.00 121.00 121.00 121.00 121.00 




ATS IDO 4 4 


..... PL A IJ C T S f S •••• rCCHNOLOOY OiVISION S K F INDUSTRIES INC. •*** P L A N C T S Y S 
rLAiicrsTs cxecuTEO usino n a s a tractiom rooel 


B E A rt 1 N T 

S Y S T E 1 

OUTPUT 

ENGLISH- 

UNITS 




STAfiE 

planet 6EAH 

SUN GEAR 

RING GEAR PLANET ORBITAL 

ROLtCII OKifTAL 


MO. 

S'»EED 

SPEED 

SPEED 

SPEED 

SPEED 

til 

f 

1 

-.lOO+OA 

.lEE^OA 

• 000 

361. 

-^591. 

€0 


STAGE 

NO. 

1 

TANGENTTAL 
GEAR LOAO* 

RADIAL GEAR TOOTH 

GEAR LOAD* TORGUE* 

Tt.9. 58A. 

IPLANCT GEAR 
INERTIA LOAD 
4A.6 


• LOAD AT ONE SEAR TOOTH. 




ROLLCII ROTATIONAL 
speco 
-. 232 » 0 ^ 
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• I* L * ?J E T S r S •••• TECHNOLOGY DIVISION S K F INDUSTRIES INC. •••• PLANET S T S ••••»• 
PLA'lCrSYS CjtECUTtO USING NASA TRACTION NOUEL 


ST3. 

1 


FATIGUE LIFEIHOURSI 

OUTER INNER single ROU OOUHLE ROU 
RACE RACE TtEARlNC BEARItlG 

•27+OS* ,2l*Qh .2S«0A .lA»OA 


H / SI SNA 

O^TER RACE 


IKNER RACE 


.05 




O O 


^ a 



tZ9 

I 

VD 


o 

c 


s 


3s 


LUlIE-LlFE FACTOR 

STj. outer RACE INNER RACE 

I .210 .210 


material factor 

OUTER RACE 

1.000 


INNER RACE 
1*000 



rENPERATURES RELEVANT TO BEARING PERFORNAHCE (OECRCCS FAHRENHEIT! 

SHAFT I. RIIIG I. RACE I« FLNG. ROLL. EL 0. FLNC. 0. RACE 0. RING BULK LUBE 
2SQ. 2S0. 250. 250. 250. 250. 250. 250. 23Cy 


r 


a. 


AT81D044 


01 -s 






p L A E T 3 Y S *•** TCCHMOLOOr DIVISION 3 K F 1NOUSTR1CS INC. pLANFfSTS ♦•••?• 

PLAHETSYS EXECUTCO USIMu NASA TRACTION NOJJEL 


B r A R J »l C S Y S T E H OUTPUT ENGLISH UNITS 

frictional heat generation RATE (BTU/llR) AND FRICTION TORDUC CLB-INl 


ST > • 

0. RACE 

I. RACE 

R.E. 

DRAG R.E.-CACC CACE>LAN0 TOTAL TORQUE 


- 

1 

'J6 .2 


Q.OOO 

3.25 0.000 A33. 10.9 

* 

• 


EIID FILH thickness* FILM REDUCTION FACTORS AND HEAT CONOUCTIVliY DATA FOR THE OUTER AND INNER RACEWAYS RESFECTIVELT 
STS. FILH INICRO ini - STARVATION FACTOR THERMAL FACTOR NEKIISCUS DIST. ilN| CONDUCTIVITY IB/HR/OEC.FI 
L .683 .507 1.00 1.00 .667 .660 5).282-02 S.A3A-82 25.3 22*8 


So 

^ ai 

s| 

S| 






AT81D044 


P L A N E T S Y S 


T^CHNOLdUY n:YlSION S K F INDUSTRIES INC 


•••• 


PLANCTSYS 


PLAiiETSYS EXECUTED IISIHG NASA TRACTION NODEL 
1) E A R I !l i SYSTCN OUTPUT ENGLISH UNITS 



LUBRICANT TEMPERATURES AND PHYSICAL PROPERTIES 


LOCATION 


DFUSITY 

VISCOSITY 

PRESSURE VISCOSITY 


<Dt!5HCES F*l 

CLn/IN3> 

KIN. riw2/seci 

DYN. itB-sec/i Nai 

COEFflCIEHT tl/PS|> 

OUTER 

2<«9.800 


•5^15*02 

•A72Q-06 

,7797-OA 

INMFR 

24<i.a00 

.3566*01 

.5^115*02 

^%720*06 

.7 797-0 A 

rfULK 

2H9.0OO 

.5366*01 

.5^15*02 

•4720-06 

,7797-OA 


CASE II A T A ENGLISH UNITS 

CATE KAIL - RING LAND. DATA CAGE SPEED DATA 

TORTJC heat ^^ATE SEP. FORCE ECCENTRICITY CAl CULATEO SPEED 
ilitC. flN-LUl (RTU/liRI IPOUHOS) RATIO IRAO/SEC) IRPMI 

I o.aoo a. 000 A.ooo-a2 o.aoo > 61.9 




• P L A N E T S Y S TECHNOLOGY DIVISION S K F INDUSTRIES INC. •••# P t A N E T S V 5 **•••• 

PLANETSYS EXECUTED USING NASA TRACTICN HOOEL 


ROLLING E L E n E N T OUTPUT FOR BEARING NUMBER i ENGLISH UNITS 


AXUMITH 
AN5LE <0E5I 

CASE 

NORMAL FORCES 
ROLLER END 

OUrrR 

INNER 

HZ STRESS 
OUTER 

15.000 

.104 

• llA 


• 031 

. .137^B5 

45.000 

.263 

.088 

.351 

• 023 

•120«05 

13.000 

.386 

• OAl 

.165 

• Oil 

•821«0A 

105.000 

• 386 

-.012 

• 000 

• 046 

.000 

135.000 

• 283 

-.059 

• 000 

.219 

.000 

165.000 

.lOA 

-.086 

• ooc 

.319 

.000 

195.000 

-.lOA 

-.086 

2A2.210 

2A2.529 

.143406 

225.000 

-•283 

-.059 

508.327 

508.547 

• 185^06 

255.000 

-.386 

-.012 

468.933 

468.979 

•179406 

255.000 

—>386 

• OAl 

A68.832 

468.678 

.179»06 

315.000 

-.283 

• 088 

502.590 

502.262 

.184406 

345.003 

-.lOA 

• llA 

221.663 

221.236 

.138406 


INNER 

t3B2tO^ 

«2G2t04 

#U?t05 
• IAUQ5 
flTG^OG 

f228«06 

•22l«0G 

•221^06 

,227^06 

«170^0G 


0. ring 
OCF t. 
•11-02 
-f39-06 

•ao-02 

.-•iO-02 

#33-0^ 

• 11-02 

•69-03 

-•11-05 

-•90-05 

-.89-03 

-• 86 - 0 % 

•91-05 


centrifugal 

FORCE 

•AA^OO 

•3A«0t 

•i6^ao 

-.A8-01 

-.23^00 

-•33^00 

-•33^00 

-•23^00 

-•AB-OI 

•IGtOO 

•3A«00 

•AA^OO 



> 

00 
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p L A H t r SYS •••• T£CHNOLOuY DIVISION S K F INDUSTRIES INC 


P L A N C T S Y S 




PLANETSYS 


A-'U’IITH 
ancle IDEGJ 

15. ana 
Aa.aOO 
7S.0J0 
liis.oon 
ijs.oaa 
1C5.000 
VJS.OUO 
."25.000 
255.00J 
p'lS.aoo 
3 15.003 
JA5.U00 


CXCCUTCO USING M A. s A TRACTION N03EL 


contact ellipse SIZE FOR STAGE NO. 1 ENGLISH UNITS 


OUTER RING 


SEMI-MAJOR 

StHI-«UN0R 

.1657-01 

•7062-03 

,.1516-01 

•6A61-03 

.1173-01 

.5025-05 

.3000 

.OuOO 

.3000 

• 0000 

.0000 

• 0000 

.13A0«00 

.6712-02 

.1715*00 

'.7313-02 

.1670*00 

.7119-02 

.16/0*00 

.7118-02 

.1703*00 

•.7285-02 

.1301*00 

.55A5-02 


INNER RING 


SEHI-HAJOR 

SEMI -MINOR 

.6714-02 

•2406-03 

.6143-02 

.2201-03 

.4778-02 

.1712-03 

.7693-02 

i2Ial.-03 

.1294-01 

.4635-03 

.1466-01 

.5253-S3 

.1338*00 

.4792-02 

.1712*00 

.6*34—02 

.1667*00 

.5971-02 

.1666*00 

.5969-02 

.1705*00 

.6109-02 

.1297*00 

.4648-02 


A 1*1 FOLLOUIHG AB ELLIPSE DIHENSION INDICATES A CONTACT LENGTH GREATER THAN ThE EFFECTIVE LENGTH 
A I*.) FOLLOVIUC AN ELLIRSE DIMENSION INDICATES A CONTACT LENGTH GREATER THAN 1.5 TIMES ThE EFFECTIVE LENGTH 
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AT81D044 


APPENDIX C 
SAMPLE OUTPUT 
♦»SKF" VERSION 


C-1 

SKF TECHNOLOGY SERViCES 

SkF INDUSTRIES. INC. . ^ 



? L A fl E T S y S TECHNOLCCY division SHF 1NOUSTRIES INC. PLAtf^TSyS ••••#• 

PLANETSYS EXECUTED USING SHF TRACTION MODEL 


THIS PLANETARY SYSTEM CONTAINS 1 STAGE. 

STAGE NO. m CONTAINS 3 IDENTICAL PLANETS* EACH SUPPORTED. BY A SPHERICAL ROLLER SEARING HAVING 2 ROH$ Of ROLLERS. 
THE MAXIMUM NO. OF RAIN LOOP ITTERATIONS ALLOUEO IS 20 AND THE RELATIVE ACCURACY REOUIREO |S .10000-07 


n 

f 


o o 

"H 71 

3 1 


O > 

30 F 

O *0 

^5 






AT81D044 




• •• 


P U A n E T S T S .... TECHMOLOCy DIVISION S K F INDUSTRIES INC. .... PLANETS TS 

PLAUETSYS EXCCUTEO USING SKF TRACTION NODCt 


O 

f 

CrI 


i 

[ 


N 0 T E S UOIHEHSIONAL U*I1TS AN5LES. .......DECREES 

DENSITY POIJNDS/CUHIC IRCH 

ELASTIC POUNDS/SQ. INCH 

FORCE. ,PQU,jQ5 

KINEHATIC Vise inches 30. /SEC. 

.........INCHES 

*'*5S POUNDS 

PRESSURE POUNOS/SQ. INCH 

SPEEDS...... ........R.P.H. 


SURFACE ROUCHNESS.............niCROtl(CHES 

temperature ....DECREES FAHRENHEIT 

THERM. -C0N0........p.T.U./lilR>F00T-0E6 F. 

21THE TERN BEARIN3 HUN3ER WHICN APPEARS BELOW IS SYNONONOUS UITH STACE NUflPCR • 

31 planetary inversion refers to THE SPECIFIC KINEMATIC CONDITION AS OESCRISEO IN THE USERS MANUAL 


> 

H 

oo 






m 
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• PL A H E T S y S ***• TCCHMOLOcr division S K F industries INC, ••*• P L A N c T s T s •••••• 

PMHETSYS CXECUrrn USING SNF TRACTIDM ROOCL 


HEAftnn ^itcm dia 

UUMSER OF flEARING 

I .2A»01 


PLAllCT 

IASS 

3.AA0 


PITCH DIA 
OF GEAR 
.395*01 


ART 


0 A I A 

V20TH OF 

planet 

1.550 


• * * • 

OIA OF PLANET NONENT PLANET TOOTH 
PL. NIRL, AXIS OF INEILUA PRESS. ANGLE 
5-2S6 ,011 2A.6 


n 

I 

4^ 


erAFtxritJ 

«UMBEH 

i 


I4U1BCR POWER THRU 

OF planets stage 

3 202.200 


plan, invers. speed of 

ISEt NOTE 31 sm 

3 1659.0 


SPEED OF 
CARRIER 
355.0 


SPECO OF 
RING 
' ^0 


pitch 01 a. 

OF SUN 

3.0 



O 


ORIGINAL PAQ2 |;| 
OF POOR QUALITY 


* P L A fl E T S Y S •••• TECHNOLOOY DIVISION S K F INDUSTRIES INC. •••• F L « N C T S T S 
PLAHETSYS executed using SKF TRICTIOH NQDEL 


hearing 

« » • « » B E A R I N C 
BEARING type 

0 A T A * < 

NO. OF 

> • • * » 

NO. OF ROLLING 

CONTjlCT 

oianctral 

nunii€H 


ROUS 

ELTS. PER ROU 

ANGLE 

clcar'ance 

1 

* SPHERICAL ROLLER BEARING 

2 

12 

15,0 

•00230 



• • • C A 0 E 0 
CAGE TYPE CAGE POCKET 

A T A ’ ♦ * • 
RAIL LAND 

rail land 

RAIL LAND 

CAGE 

•IU1BCR 

clearance 

WIDTH 

diaheter 

CLEARANCE 

HASS 

l 

ROLLING ELEMENT HTUING .0100 

• 0000 

• 000 

• 000 

• 01 


nrARjur. 

iiuMuni 

1 


STEEL TYPE* 
INNER RING 
1-SO FOR-7 


• steel 0 a T a • » • 

LIFE steel type* 

FACTOR OUTER RING 

1.000 CARBURIZED 9310 


LIFE 

FACTOR 

l.OOO 
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E . 1 BACKGROUND 

The temperature portion of the Program is designed to 
produce temperature maps for an axisymmetric mechanical system 
of any geometrical shape. The mechanical system is first ap*" 
proximated by an equivalent system comprising a number of ele- 
ments of simple geometries. Each element is then represented 
by a node point having either a knovm or an unknown temperature. 
The environment surrounding the system is also represented by 
one or more nodes. With the node points properly selected, the 
heat balance equations can be set up accordingly for the nodes 
of unknown temperature. These equations become non-linear when 
there is free convection and/or radiation between two or more of 
the node points considered. The problem is therefore reduced to 
solving a set of linear and/or non-linear equations for the same 
number of unknown nodal temperatures. It is obvious that the 
success of the approach depends largely on the physical subdi- 
vision of the system* If the subdivision is too fine, there 
will be a large number of equations to be solved;, on the other 
hand, if the subdivision is too crude, the results may not be 
reliable. 

In a system consisting of rolling bearings, for the sake of 
simplicity, the elements considered are usually axially symmet- 
rical, e.g., each of the bearing rings can be taken as an ele- 
ment of uniform temperature. For an element which is not axially 
symmetrical, its temperature is also assumed to be uniform and 
its presence is. assumed not to distort the uniformity in tem- 
perature of a neighboring element which is axially symmetrical. 
That is, the non- symmetrical element is represented by an 
equivalent axially symmetrical element with approximately the 
same surface area and material volume. This kind of approxima- 
tion may seem to be somewhat unrealistic, but with properly de- 
vised equivalent systems, it can be used to solve complicated 
problems with results satisfying most of the important engineer- 
ing requirements . 

The computer program can solve the heat-balance equations 
for either the steady state or the transient state conditions 
and produce temperature maps for the mechanical system when the 
input data are properly prepared. 

E.2 BASIC EQUATIONS 

E.2.1 Heat Conduction 

The rat© of heat flow qci,j(W) that is conducted from node 
i to node j may be represented by. 
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^ij 


- 




ti and tj are the temperatures at i and j, respectively, Aj,, j 
the area"^ normal to the heat flow (m^) , Lm the distance (mf and 
the thermal conductivity between i and j, (W/M®C) . 

Assuming that the structure between point i and j is 
composed of different materials, on equivalent: heat conductivity 
may be calculated as follows: 



Pig. E-1 


^ 1^1 ^2 ^2 


^ij ” ^1 ^2 



E.2.2 CONVECTION 

The rate of heat flow that is transferred between a solid 
structure and air by free convection may be expressed by 


E“3 
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^yi,j • *i,j 


where 


SIGN- 1 

■ Ulr if (ti* - tj). < 


0 

t j ). < 0 


in which 


2.5 • 10 ^ W/m^- (degC)^*^^ for hot surfaces facing upward 


and cold surfaces facing downward 

®ij"^1.4' 10"^ W/m^- (degC)^*^^ for hot surfaces facing downward 

and cold surfaces facing upward 


tl.8‘ 10“^ W/m^- (degC)^*^^ vertical surfaces 


For other special conditions , ^ ij must be estimat.ed hy referring 
to heat transfer literature. 


The rate of heat flow that is transferred between a solid 
structure and a fluid by forced convection my be expressed by 

^i,j " ^ 


in which is the heat transfer coefficient 


Now, with a ■a- • introduce the Nusselt number 

*J » 


the Reynolds number 


N - ^ 

A 


and the Prandtl number 


K ■ UL 
® “u 


p a P Cp 

where 

L is a characteristic length which is equal to the diameter 
in the case of a cylindrical surface and is equal to the 
plate length in case of a flat surface (m) 
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U is a characteristic velocity which is equal to the difference 
between the fluid velocity at some distance from the surface 
and. the surface velocity (m/sec) 

X is the fluid thermal conductivity (W/M®C) 

2 

\) is the fluid kinematic viscosity (M /sec) 
p is the fluid density (kg/m^) 

C is the fluid specific heat (J/kg°C) 

P 

For given values of R. and Pf the Nusselt number Nu and 
thus the heat transfer coefficent may be estimated from one 
of the following expressions; 

Laminar flow along a flat place; R, <2300 

- 0.323 

Laminar flow of a liquid in a pipe: 

Na ■ K ■ 

where D is the pipe diameter and L the pipe length 
Turbulent flow of a liquid in a pipe: 

- 0.027 . Rg^’®. 

Gas flow inside and outside a tube; 

,57 

Liquid flow outside a tube : 

Nu ■ 

Forced free convection from the outer surface of a 
rotating shaft 

\ » 0.11 Co . 5 Rg ^ 

O 

¥&?i?iStl® J^ey^olds number e is developed by the shaxt 

» M 7T 
V 
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ia which o) is the angular velocity (rad /sec) 

D is the shaft diameter (m) 

The ‘average coefficient of forced convection to the lubri- 
cating oil within a rolling contact bearing may be approximated 
by, 

k 

am 0.0986 |~C1 + 

using for outer ring rotation 
• for inner ring rotation 

in which N is the bearing operating speed (rad/sec) 

D is the diameter of the rolling elements (mm) 

“m is the bearing pitch diameter (mm) 

B is the bearing contact angle (degrees) 


E.2.3 fluid flow 


The rate of heat flow that is transferred from fluid 
node i to fluid node j by fluid flow is 


‘»£i. j • 1’'^ 


13 


S 


V. • is the volume rate of flow from i to j , It must be observed 

that the continuity of mass requires the following equation to 
be satisfied 


provided the fluid density is constant. The summation should 
be extended over all nodes i within the fluid which have heat 
exchange with node j by fluid flow. ^ 

E.2.4 HEAT RADIATION 

The rate of heat flow that is radiated to node j from node 
i is expressed by 

where 

tj * Temperature of node j in °C 
= Temperature of node i in °G 
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and the value of the coefficient ^i,j depends on the geometry 
and the emissivity or the absorptivity of the bodies involved. 

For radiation between large, parallel and adjacent surfaces 
of equal area, ^i,j and emissivityi c l',j,6 i,j is obtained from 
the equation 

*i,J ®*i,J 

where T , the Stefan ’•Boltzmann Constant, is 
0 . 5.76 • 10"® W/m^ /CdegK)^ 

For radiation between concentric spheres and coaxial cylind»’ 
ers of equal emissivity, ®i,j, ®i, j is given by the equation 



®iij ^^i.i 


A. * 

®i 


1 + ( 1 - 


A 

where a is as above ^ i,j is the area of the enclosed body 
and A*i,j is the area of the surrounding body, i.e.A- i <A'\. .. 

A I J ^ > J 

Expressions for <Si,j that are valid for more complicated 
geometries or 'for different emissivities may be found in the 
heat transfer literature. 


E.2.5 CALCULATION OF AREAS 

In the case of conduction heat transfer in the 
direction A^ ^ is given by the equation (Pig. E-3) 

^m’ 


Referring to the input instructions, Section 2.4.R 

+ 

T 


1 m 


ri + rj 


Ar * 
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In the case of heat transfer in the radial direction, i 
. is obtained from the expression 

H,) ■ ^1 ' V 4 ■ « 

: f 

and similarly for the radiation term above (Figure B-4 (c)) 

A^ 4 4 ■ 2ff r*^H 
i»J ™ 


4 - H 

in which H is the length of the cylindrical surface! where 'heat 
is conducted between i and j, r is given by the same equation 
as above (Fig. E-4 (a) ); where"^heat is convected between i and 
r is the radius of the cylindrical surface CFig. E ".4 (b) ); 
wRere heat is radiated between i‘ and j, r is the radius of the 
enclosed cylindrical surface and r ^ the radius of the surround- 
ing cylindrical surface (Fig. E-4 ic) ) 
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E.3.1. transient analysis 

For the transient analysis all of the data pertaining to 
the node to node heat transfer coefficients must be provided by 
the input. Additionally, the volume and the specific heat at 
each node is required. 


E-9 

SKF TECHNOLOGY SERVICES 

SKF INDUSTRieS, :NC. 



AT81D044 


PART 2; EQUATIONS USED IN TEMPERATURE CALCULATIONS 
E". 4 Temperature Calculations 


Subsequent to each calculation Qf bearing generated heat 
rates, either the steady state or transient temperature mapping 
solution scheme may be executed. This set o£ sequential calcu- 
lations is terminated as follows : 

1. For the steady^state case; when eaqh system temperature 
is within EPl ^Centrigrade of its previously pre- 
dicted value, EPi is specified by the user* If it is 
zero or left blank, a default value of 1° Centigrade is 
used. This criteria implies that the Steady state 
equilibrium conditions had been reached, 

2. The transient calculation terminates when the user' 
specified time up is reached or when one of the system 
temperatures exceeds 600*^C. 

E , 4 . 1 . Steady State Temperature Map 

The mechanical structure to be analyzed is thought of as 
divided into a number of elements or nodes, each represented by 
a temperature. The net heat flow to node i from the surrounding 
nodes j, plus the heat generated at node i, must numerically 
equal zero. This is true for each node i, i going from 1 to n, 
n being the number of unknown temperatures. 

After each calculation of bearing generated heat, which 
results from a solution of the planet bearing system portion 
of the program, a set of system temperatures is determined 
which satisfy the system of equations; 

q^ • q^^ q^^ ■ 0 for all temperature nodes i (Bl) 

W K G 1* S Q 

^oi is the heat flow from all neighboring nodes to node i 

*^gi is the heat generated at node i. These values may 

be input or calculated by the planet bearing program 
as bearing frictional heat 
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This scheme is solved with a modified Newton- Raphson 
method which successful!/ terminates when either of two 
conditions are met: 



< E)P2 for all nodes i 


CB.i) 


lyhere: At represents the Newton-Raphson correction to the 

temperature t at a given iteration such that, 

^N+1 " ■► At and N •► 1, and N, refer to the 

next and current iteration respectively. 


EP2 is a user specified constant. If EP2 is left blank 
or set to zero (0) a default value of 0.001 is used. 


A second convergence criterion dependent upon EP2 is also 
used. In the system of equations, q . •► q . » 0 for all nodes 
i, absolute convergence would be otalned it^the right hand 
side (EQ) in fact reduced to zero (0). Usually a small residue 
remains at each node, such that (q^^^ * CEQ)i. 

The second convergence criterion is satisfied if 


n 

i=l 


CEQ) 


n 


1 100 X EP2 


(E.3) 


where n ■ number of equations in thermal solution 
= number of unknown temperatures 
B.4.2 Transient Temperatures 


In the transient case the net heat q. transferred to a 
node i heats the element. It is thus necessary for heat balance 
at node i that the following equations are satisfied. 


Pi 


Cp. 




dti 

3TT 


“ q. 


(E.4) 
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where P 

• 

» density 

s 

■ specific heat 

V 

■ volume of the element 

■t; 

T 

■ temperature 
• time 


The temperatures, t^- at the time of initiation T ■ T are 
assumed to be known, ^ that is ^ 

t. (T.) • t . i» 1,2,».,, n (B..S) 

9 OX 


The problem of calculating the transient temperature f 

distribution in. a bearing arrangement thus becomes a problem | 

of solving a system of non-linear differential equations of ; 

the first order with certain initial values given. The equations 
are non-linear since they contain terms of radiation and free I 

convection, which are non-linear with temperature as will .be 
shown later. The simplest and most economical way of solving 
these equations is to calculate the rate of temperature in 
crease at the time T * Ti. from equation e* 4 and tj|ien calcu- 
late the temperatures at*^time Tj^ + AT from 

*k+l • *k * — * ‘k * ‘’it 4T CE-6) 

dT pCpV 

If the time step <xT used as program input is chosen 
too large, the temperatures will oscillate, and if it is chosen 
too small the calculation will be costly. It is therefore 
desirable to choose the largest possible time step that does 
not give an oscillating solution. The program optionally cal- 
culates such a time step. The step is obtained from the condi- 
tion, 

^^i. k»l 

dt^ ^0 i - 1 , 2 , . . . , n (E . 7 ) 

If this derivation were negative, the implication would 
be that the local temperature at node i has a negative effect 
on its future value . This v/ould be tantamount to assefting 
that the hotter a region is now, the colder it will be after 
an equal time interval. An oscillating solution would result. 
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Differentiating equation (E. 6) for node i, one has as 
condition. 


‘^^i.k-*-! - 1 + ^'^i . ^ 

dt. ^ ^ iSi'^i ‘^^i, k 

XplS. 


i ■ 1, 2,...n (E. 8) 


The derivative dq^^ k/*^^i k calculated numerically 

‘’l.k'H - ‘'i,k 

For each node the value of aT. giving a value of zero to, 
the right hand side of Eqn. (E.8)^is calculated.' 


A value of aT rounded off to one significant digit smaller 
•than the smallest of the aT. given by Eq. OB. 8) is used. 

If the transient thermal scheme is being used interactively 
with the bearing subprograms the user must specify a small enough 
time step between calls to the bearing subprograms in order that 
the variation in bearing generated heats, with time, accurately 
reflects the physical situation. At first a trial and error 
procedure will be required to effectively use the program in this 
mode, however, experience will increase the user’s effectiveness. 

E. 4.3 Calculation of Heat Transfer Rate 


The transfer of heat within a medium or between two media 
can occur by conduction, convection radiation and fluid flow. 

All these types of heat transfer occur in a bearing appli 
cation as the following examples show*. 

1. Heat is transferred by conduction between inner ring 
and post and between outer ring and housing. 

2. Heat is transferred by convection between the surface 
of the housing and the suTrounding air. 
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3. When the bearing is lubricated and cooled by circulating 
oil, heat is transferred by fluid flow. 

Therefore, in calculating the net flow to a node all the 
above mentioned modes of heat transfer will be considered. 


E4,3.1. Generated Heat 

There may be a heat source at node i giving rise to a heat 
flow, to be added to the heat flowing from the neighboring nodes. 


In the case that the heat source is a bearing, it may 
either be considered to produce known amounts of power, in which 
case constant numbers are entered as input to the program, or 
the planet bearing program ma.)^ be used to calculate the bearing 
generated heat as a function of bearing temperatures. 


E . 4 . 3 . 2 . Conductio n 

— M i - . 1 w»ii^M= =g=g= t 

The heat flow q . • which is transferred by conduction 
from node i to node j,'is proportional to the difference in 
temperature (t- - t.) and the cross-sectional area A and 
is inversely proportional to the distance 1 between the two 
points, thus 


lei, j “ ^^i " 


CE.10) 


where ^ » the thermal conductivity of the medium. 
E. 4. 3.3 Free Convection 


Between a solid medium such as a metallic body and a 
liquid or gas, heat transfer is by free or forced convection. 
Heat transfer by free convection is caused by the setting in 
motion of the liquid or gas as a result of a change in density 
arising from a temperature differential in the medium. With 
free convection between a solid medium and air, the heat 
energy q . . transferred between nodes i and j can be calculated 
from the' equation, 

’’vi.j " l*i' . CE.H) 

where the film coefficient of heat transfer by free 
convection 

A = the surface area of contact between the media 
d » is an exponent, usually = 1.25, but any value 
can be specified as input to the Program 


SKF TECHNOLOGY SERVCES 

SKF INDUSTRIES. 


E-14 


ORIGINAL PAGE IS 
OF POOR QUALITY 


I; e 


AT8ID04'4 


SIGN = 


. Tl if 

if t^<tj 


The last factor is included to give the expression q„. . 
a correct sign. 


The value of a can be calculated for various cases, see 
syacob and Hawkins, {16}. 

B.4,3.4. Forced Convection 


Heat transfer by forced convection takes place when liquid 
or gas moves around a solid body, for example, when the liquid 
is forced to flow by means of a pump or when the solid body is 
moved through the liquid or gas. The heat flow q^i j transferred 
by forced convection can be obtained from the follo^ng equation; 

^wi,j * (E»12) 

* 

where ^ is the film coefficient of heat trasfer during 
forced convection. This value is deisendent on 
the actual shape, the surface condition of the 
body, the difference in speed, as well as the 
properties of the liquid or gas. 

In most cases, it is possible to calculate the coefficient 
of forced convection from a general relationship of the form, 

(E.13) 


where a, b, and c are constants obtained from handbooks, 

such as (17} , R and P are dimensionless numbers 
defined by e r 

Nu ■ Nusselt's number = owL/x 
L ■ characteristic length 
X ■ conductivity of the fluid 

R ■ Reynold's number * ULp/n 

U * characteristic speed 

p » density of the fluid 

n * dynamic viscosity of the fluid 

P^» Prandtl's number = 

Cp* specific heat 
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The program caji Use a constant value of the coefficient of conyec- 
tioHi or let it vary with actual temperatures, the variation 
being determined by how the viscosity varies. Input can then be 
given in one of three ways, for each coefficient. 

Constant viscosity 


1* Values of the parameters o£ equation rB.lS) 

are given as input and a constant value of 
is calculated by the program. 

Temperature dependent viscosity 

2. The coefficient for turbulent flow and 
heating of petroleum oils is given by 

®w " ^9 * n (t) ^10 

^Vhere; kg and kj_Q are given as input together 

with viscosity at two different temperatures. 

3. Values of the parameters of equation (e» 13) 
are given as input. Viscosity is given at 
two different temperatures. 


E.4.3.S Radiation 

If two flat parallel, similar surfaces are placed 
close together and have the same surface area A, the heat 
energy transferrc?d by radiation bet\^een nodes i and j 
representing those bodies, will be, 

Where: e is the surface emissivity. The value of 

the coefficient s is an input variable and 
varies between l for a completely black surface 
and 0 for an absolutely cle|.n surface. In 
addition cr is Stefan-Bolttmann’ s radiation 
constant which has the value 5.76 x 10 “® watts/ 
m2 - (*K)4 and t- and tj are the temperatures ,Oj^n 
at points i and j . • 

Heat transfer by radiation under othex' conditions can 
also be calculated, {16}. The following equation, for 
instance, applies between two concentric cylindrical surfaces 
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q -e g *i rfh * 273)* - f^1 * 271^^1 

1. > (1-e) (A^/A,) 


(P .16) 


Where A,^ is the aree o£ the inner cylindrical surface 
is the area of the outer cylindrical surface 

■C* 

a 4, 3.6 Fluid Flow 


Between nodes established in fluids, heat is transferred 
by transport of the fluid itself and the hciat it contains. 
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Figure E-5 Convective Heat Transfer 

Figure E -S shows nodes i and j at the midpoints of consecu 
tive segments established in a stream of flowing fluid. 

The heat flow j through the boundary between nodes 

i, and j can be calculated as the sum of the heat flow 
through the middle of the element i, and half the heat 
flow q^,^ transferred to node i by other means, such as 
convection. 

The heat carried by mass flow is, 


■Ifi ■J>i Si ’'i ■ ’'i 


Where Vj|_ ■ the volume flow rate through node i 


CE.17) 


The heat input to node i is the sum of the heat generated 
at node 1 (if any) and the sum over all other nodes of the 
heat transferred to node i by conduction, radiation, free 
and forced convection. 
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‘^ol " ■*' ^ ‘^vij ^ ^wi,j 


(B.18) 


The heat flow between the nodes of Fig. E*S is then» 


V.J ’ ’«i * 


(E.19) 


If the flow ftom node i is dividing between nodes j and k, Figure 
E« 6 then the h^Kat flow is calculated from 


■»Ui,j ■ ■'il (<’fi ♦ 


(B.20) 


IVhere ■ the proportion of the flow at i going to 

node j, O^K^j i 1* i* specified at input 


K.-K. .K. 

] ij 1 

Ki.* k»i,* ki 


K.» K.+Ku 
1 3 K 



FIGURE E -6 

DIVIDED FLUID FLOW FROM NODE i 
E.4.3.7 Total Heat Transferred 


The net heat flow rate to node i can be expressed as 


IIA » 

'll - '’G.i * * ^ui.j * “'vi.i ♦ * ^Rij) 
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The summation should include all nodes j , both with unknown 
temperatures as well as boundary nodes, at which the 
temperature is known, so long as they have a direct heat exchange 
with node i. 


This expression is a non-linear function of temperatures 
because of the terms qw qR . Therefore the equations to 
be solved for a steady state solution are non-linear* The 
subprogram SOLVXX for solving non-linear simultaneous equations 
is used for this purpose . 

E'.4.4 Conduction Through a Bearing 

As described in Section E.4.3.2 the conduction between two 
nodes is governed by the thermal conductivity parameter X of 
the medium through which conduction takes place. The value ’ 
of X is specified at input. 

An exception is when one of the nodes represents a bearing 
ring and the other a 'set of rolling elements. In this case 
the conduction is separately calculated using the principles 
described below. • 

E . 4 . 4 . 1 Thermal Resistance 

It is assumed that the rolling speeds of the rolling elements 
are so high that the bulk temperature of the rolling elements 
are the same at both the inner and outer races, except in a 
volume close to the surface. The resistance to heat flow can 
then be calculated as the sum of the resistance across the 
surface and the resistance of the material close to the 
surface . 


The resistance is defined implicitly by 

At ■ n *q (E.22) 

where 

At is temperature difference 
q is heat flow 

The resistance due to conduction through the EHD film is 
calculated as 

a - h/ CXA) (E .23) 
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where h is taken to be the calculated plateau film 
thickness 

A is the Hertzian contact area at the specific 
rolling element^ring contact under considera* 
tion 

X is the conductivity of the oil. 

The geometry is shown in Figure ' e- 7 (a). 


So far, a constant temperature difference between the 
surface has been assumed. But during the time period of 
contact, the difference will decrease because of the finite 
thermal diffusivity of the material near the surface, Fig. E -7(b). 

To points at a distance from the surface this phenomenon 
will have the same effect as an additional resistance O 5 acting 
in series with 0^ . 

This resistance was estimated in {ig} 


« 1 t 

% -Xl ^ 2b .V 

re, i 

wherejl^Q -contact length, or in the case of an 

elliptical contact area, 0.8 times the 
major axis 

X - heat conductivity 
f - thermal diffusivity X/ (pCp) 
p - density 
Cp - specific heat 



b - half the contact width 
V - rolling speed 
The resultant resistance is 
^res ■ ^1 + ^2 


E-21 


1 SKF TECHMDirS/ SERVICES 

SkF INDUS"” 


CE.24) 


(E. 25) 


ORIQINM. PAOe n ATaiD044 

OF POOR OUAUTf 


There is one such resistance at each rolling element. 
They all act in parallel. The equivalent resistance eqv. 
is thus obtained from 


1 mi 1 ^■; 26 ) 

^ ©qv ^ ^‘'^resRi 

I 

E..4.S Consideration of Multiple, Axially Symmetric 

Subsystems within the Primary Thermal System. 

The temperature mapping heat dissipation analysis is built 
on the assumption that the system is basically axially symmetric- 
al so that each temperature node is essentially a circular ring 
with no temperature deviations around the circumference. 

Solutions of the planetary system problem assumes that each 
nf the planets within a stage transmits equal power and is 
subjected to identical load and speed conditions. This assump- 
tion permits the calculation of frictional heat for one planet 
bearing (one or two rows of rollers) but allows us to 
attribute the same frictional heat to each of the planets 
within the stage. The heat flow within a given planet is 
governed by its axially symmetric heat transfer characteristics. 

The heat flow in the primary system is affected by the heat 
generated by all individual planets and its symmetry about 
another axis . To account for this the heat transfer coefficients 
between nodes in the prime and subsystem are direct and so that 
when heat flows from thsprimary system to the subsystem only 
(1/Npl) of the heat enters the subsystem where Npl is the number 
of planets. The inverse is true as heat leaves a subsystem and 
enters the primary system. 

The subsystem-primary system heat transfer interaction is 
taken into account with program input, on card type T8 . 
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P-1 . INTRODUCTION 

This Appendix is a supplement to the User's Manual for SKF 
Computer Program PLANETSYS. It highlights the differences between 
the PLANETSYS/SKF and PLANETSYS/NASA program versions. In Section 
F-2 below, the differences between the mathematical models used in 
the two versions of the program are stated explicitly. These dif- 
ferences are related to: 

1) The calculation of EHD film thickness 

2) Concentrated contact traction 

In section F-3, the differences between the two versions in terms 
of program execution are specified. 

F-2. MATHEMATICAL MODELS 


F-2.1 EHD Film Thickness 


PLANETSYS/SKF uses the Dowson Higginson (13) equation to 
calculate film thickness for line contact. The calculated value 
is multiplicatively modified by a film thickness thermal reduc- 
tion factor developed by Cheng (8) and by a starvation reduction 
factor developed by Chiu (9) . 

PLANETSYS/NASA uses the equation developed by Loewenthal 
et al (21) in calculating film thickness, for both point and line 
contact. 

F-2. 2 Concentrated Contract Traction 

The concentrated contact traction model used in PLANETSYS/ 
SKF accounts for lubricant shear and asperity interaction. A 
semi- empirical model (10), is used to calculate an EHD lubricant 
shear coefficient. Asperity effects are introduced by determin- 
ing the portion of the contact load carried by the asperities, 
using the analysis of Tallian (11), and then calculating the re- 
sulting traction as the product of the normal load carried by the 
asperities times the asperity friction coefficient. In equation 
form the traction force is: 


p ^ 2 
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!eHD ’^EHD ■*" .^ASP ^ASP 

(F-l) 

!eHD ^ *^ASP 

(F-2) 


Q is the normal load 

QgHD the normal load carried by the EHD film 

lipun friction coefficient which develops from 

lubricant shear 

Qa 3 p is the normal load carried by asperities 
V 1 a 3 p is the asperity friction coefficient 
F is the traction force 


PLANETSYS/NASA calculates concentrated contact traction across 
the EHD film only, according to the model developed by Allen, et 
al (22)* This model calculates the traction force by first calcu- 
lating the shear stress according the the Newtonian fluid shear 
equation* 


T 



(F-3) 


where x is the shear stress 

n is the dynamic viscosity 
V is the surface relative sliding velocity 
h is the lubricant film thickness 

The lubricant viscosity is assumed to be an exponential func 
tion of pressure of the form: 


n 



(F-4) 


where is the dynamic viscosity at atmospheric pressure 

a is the pressure viscosity coefficient defined 
implicitly above 

P is the normal pressure 
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Allen requires that the shear stress not exceed a specified frac 
tion of the normal stress such that 


T - n J if n ^ < r, 

T ■ fp if n ^ > fP and n ^ 

where, is the critical shear stress for which a value 
of 0.0069 N/mm^ (1000 psi) is used. 

f is called the lubricant friction coefficient and has 
been determined for specific lubricants. Typical values of 

f lie in the range 0.05 f < 0. 08. 

• 

Having calculated the shear stress, the traction force is 
obtained by integrating the shear stress over the respective con- 
tact area. 


(F-S) 
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F-3. SKF VS. NASA VERSIONS - PROGRAM USE 

The selection of the desired PLANETSYS version has been made 
possible by the inclusion of two separate Map Statements for the 
Univac 1100 computer: 

1) The original map statement 
0MAP,S R,A 

enables an execution with the SKF film thickness and 
traction models, Subroutines FMIXPL, EHDSKF, and 
FRICTN are called to compute the EHD traction forces. 

2) A new map statement 
0MAP.S RNASA,ANASA 

enables the use of the NASA film thickness 
and traction models, Subroutines FMIXPL/NASA, 

FILM, and ALLEN are called to compute the EHD 
traction forces. 
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